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A theoretical framework for low-frequency electromagnetic (drift-)kinetic turbulence in 
a collisionless, multi-species plasma is presented. The result generalises reduced magne¬ 
tohydrodynamics (RMHD) and kinetic RMHD (Schekochihin et al. 2009) to the case 
where the mean distribution function of the plasma is pressure-anisotropic and different 
ion species are allowed to drift with respect to each other—a situation routinely encoun¬ 
tered in the solar wind and presumably ubiquitous in hot dilute astrophysical plasmas 
such as the intracluster medium. Two main objectives are achieved. First, in a non- 
Maxwellian plasma, the relationships between fluctuating fields (e.g., the Alfven ratio) 
are order-unity modified compared to the more commonly considered Maxwellian case, 
and so a quantitative theory is developed to support quantitative measurements now 
possible in the solar wind. Beyond these order-unity corrections, the main physical fea¬ 
ture of low-frequency plasma turbulence survives the generalisation to non-Maxwellian 
distributions: Alfvenic and compressive fluctuations are energetically decoupled, with the 
latter passively advected by the former; the Alfvenic cascade is fluid, satisfying RMHD 
equations (with the Alfven speed modified by pressure anisotropy and species drifts), 
whereas the compressive cascade is kinetic and subject to collisionless damping (and for 
a bi-Maxwellian plasma splits into three independent collisionless cascades). Secondly, 
the organising principle of this turbulence is elucidated in the form of a conservation law 
for the appropriately generalised kinetic free energy. It is shown that non-Maxwellian 
features in the distribution function reduce the rate of phase mixing and the efficacy 
of magnetic stresses, and that these changes influence the partitioning of free energy 
amongst the various cascade channels. As the firehose or mirror instability thresholds 
are approached, the dynamics of the plasma are modified so as to reduce the energetic 
cost of bending magnetic-field lines or of compressing/rarefying them. Finally, it is shown 
that this theory can be derived as a long-wavelength limit of non-Maxwellian slab gy- 
rokinetics. 
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1. Introduction 

Reduced magnetohydrodynamics (RMHD) is a non-linear system of fluid equations 
used to describe anisotropic fluctuations in magnetised plasmas at lengthscales L much 
larger than the ion gyroradius pt and at frequencies w much smaller than the ion gyrofre- 
quency It was initially used to model elongated structures in tokamaks (Kadomtsev 
& Pogutse 1974; Strauss 1976, 1977) but has since become a standard paradigm for 
astrophysical contexts such as solar-wind turbulence (Zank & Matthaeus 1992a, 6; Bhat- 
tacharjee et al. 1998) and the solar corona (Oughton et al. 2003; Perez & Chandran 
2013). 

Although RMHD was initially derived from incompressible ideal MHD, a collisional 
fluid theory, it can also be obtained without assuming the plasma to be collisional 
(Schekochihin et al. 2009, hereafter S09). The resulting set of fluid-kinetic equations 
describing both Alfvenic (RMHD) and compressive, i.e., density and magnetic-field- 
strength, fluctuations is referred to as kinetic reduced magnetohydrodynamics (KRMHD). 
S09 argued that KRMHD is an appropriate description for small-scale solar-wind fluc¬ 
tuations, which are anisotropic (e.g., Horbury et al. 2012) and weakly collisional (e.g., 
Bruno & Carbone 2005) , as well as for inertial-range turbulence in the hot ionised phase 
of the interstellar medium and in the intracluster medium of galaxy clusters. 

Two of the assumptions of KRMHD are that the equilibrium distribution functions 
of all species are Maxwellian (and, therefore, that the equilibrium pressure is isotropic) 
and that there is only one ionic species. The former assumption works well for plasmas 
such as the interstellar medium, where collisions are weak (A m f p pt and va -C fli, 
where A m f p is the collisional mean free path and Va is the ion-ion collision frequency) 
but non-negligible (A m f p <C L and Va w). However, the collisional mean free path in 
space plasmas is on the order of 1 au —the distance between the Sun and the Earth—and 
proton (H ions), alpha (He ions), and electron pressures in the solar wind are observed 
to be highly anisotropic with respect to the local magnetic-field direction (e.g., Hellinger 
et al. 2006; Stverak et al. 2008; Bale et al. 2009; Maruca et al. 2012). The observed 
distribution functions in the solar wind (especially the electron one) also exhibit non- 
Maxwellian suprathermal tails (see Maksimovic et al. 2005; Marsch 2006, and references 
therein) containing small (~5% of the total density) populations of energetic particles. 
In the intracluster medium, where A m f p ~ 0.1-30 kpc is many orders of magnitude larger 
than pi ~ 1 npc, conservation of particles’ first adiabatic invariant during (macroscale) 
turbulent stretching of the magnetic field is expected to render the distribution function 
anisotropic (e.g., Schekochihin et al. 2005; Schekochihin & Cowley 2006; Kunz et al. 
2011). How such anisotropic distribution functions affect the turbulent cascade in these 
systems is presently unknown, and it seems dangerous to describe their dynamics with 
a set of equations built upon the assumption of isotropy. 

The assumption of a single ionic species is equally unwarranted in the solar wind, where 
the abundances of alpha particles and heavy ions have been established observationally 
for nearly fifty years (for a review, see von Steiger et al. 1997). The protons and alphas 
(as well as many other ions) drift with respect to the centre-of-mass frame, often with 
parallel drift velocities on the order of the Alfven speed (Asbridge et al. 1976; Marsch 
et al. 1981, 1982; Neugebauer et al. 1994; Goldstein et al. 2000). 

In this Paper, we present a generalisation of KRMHD to account for non-Maxwellian 
distribution functions and multiple ionic species. Our theory of pressure-anisotropic 
KRMHD can be applied to a broad range of plasmas that exhibit these characteris¬ 
tics and satisfy the original assumptions of KRMHD: that the turbulent fluctuations 
are small compared to the mean field, are spatially anisotropic with respect to it, have 
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frequencies small compared to the ion cyclotron frequency, and have lengthscales large 
compared to the ion gyroradius. The purpose of this Paper is twofold: first, to explain 
what is the organising principle governing the turbulent cascade in a pressure-anisotropic 
plasma ; and second, to provide a quantitatively correct set of equations that describe that 
turbulence. A reader interested only in the former can proceed directly to Section 5.1, 
where we present a generalised free-energy invariant describing a turbulent cascade from 
large to small scales (as well as into phase space) in a non-Maxwellian plasma. A reader 
interested only in the latter will find those equations summarised in Section 2.6. 

Before proceeding with the derivation, we caution that pressure-anisotropic plasmas are 
subject to a variety of kinetic microscale instabilities if their pressure anisotropy p± — p\\ 
becomes larger than the magnetic pressure (times a factor of order unity). Some of these 
instabilities (e.g., ion-cyclotron, whistler) are ordered out of KRMHD by its restriction 
to sub-Larmor frequencies. Others, namely firehose and mirror, are included (§§3.1, 4.4), 
but have growth rates that increase without bound with wavenumber due to the exclu¬ 
sion of finite-Larmor-radius effects that would have regularised them. Plasmas exhibiting 
super-Alfvenic inter-species drifts can also be subject to cyclotron- and Landau-resonant 
electromagnetic instabilities (e.g., Marsch & Livi 1987; Gary 1991; Daughton & Gary 
1998; Verscharen et al. 2013). With these complications borne in mind, KRMHD as a 
quantitative theory is only suitable for kinetic turbulence residing within the microscale 
stability boundaries. We will show that pressure anisotropies and interspecies drifts, even 
those lying within the stability boundaries, lead to order-unity modifications of the rela¬ 
tions between different fluctuating fields and of the kinetics of the phase-space turbulent 
cascade compared to what was deduced previously for a Maxwellian plasma. 

The Paper is organised as follows. In Section 2 we derive the generalised KRMHD equa¬ 
tions starting from the drift kinetics. Several consequences of these equations are detailed 
both mathematically and physically there and in Sections 3 and 4. These include linear 
waves and stability (§§3.1, 4.4, Appendix B), nonlinearly conserved quantities (§§3.3, 
4.1, 5.1), and their combined implications for the efficacy of phase mixing (§§4.5, 4.6). In 
Appendix C, we re-derive KRMHD systematically from a generalisation of slab gyroki- 
netics (Howes et al. 2006; S09) to non-Maxwellian distribution functions and multiple 
ionic species. The latter approach enables analytical and numerical studies of fluctua¬ 
tions at and below the ion gyroradius, including the effects of pressure anisotropy on 
the nonlinear perpendicular phase mixing and the phase-space cascade of kinetic Alfven 
waves and entropy. These topics will be the subject of a separate publication. Finally, in 
Appendix E we provide a list of frequently used symbols and their definitions. 

2. General Non-Linear Equations of KRMHD 

2.1. Kinetic MHD 

We begin with the equations of kinetic MHD (KMHD), as derived by Kulsrud (1964, 
1983). KMHD is a hybrid fluid-kinetic theory, appropriate for scales A; -1 larger than the 
ion gyroradius (kpi <C 1) and frequencies w smaller than the ion gyrofrequency (w < flj) 
in a weakly collisional (pi -C A m f p ) plasma. The equations for species s are the continuity 
equation 

n s = -n s V-u s , (2.1) 

the momentum equation 

m s n s ^ + u s ■ u s = -V • P s + q s n s (-E + ^ x B^j + F s , (2.2) 


d_ 

dt 
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and the drift-kinetic equation 

D/s D In B w± df s / g s E\\ w\ - D u ±s A df s = / d/ a \ 

Df D t 2 \ m s 2 D t ) dv\\ \ dt J c 


(2.3) 


where D/D t = d/dt + u± s • V + t?|| b • V. Our notation is standard: m s , q s ( =Z s e ), n s , u s , 
and P s are, respectively, the mass, charge, number density, mean velocity, and pressure 
tensor of species s; B is the magnetic field; and E is the electric field. The rate of change 
in the momentum of species s due to inter-species collisions, F a , is obtained by taking the 
first velocity moment of the collision term (df s /dt) c . The distribution function of species 
s, f s = fa(v\\,w±), is written in terms of the velocity-space variables dm and w± , which 
are measured parallel and perpendicular, respectively, to the magnetic field direction, 
b = B/B. As in Kulsrud’s original formulation, we use as the perpendicular kinetic 
variable — u± s , the perpendicular velocity peculiar to the mean perpendicular 

flow of species s. The parallel component of Equation (2.2) is redundant, as it may be 
straightforwardly obtained by taking the first parallel velocity moment of the drift-kinetic 
equation (2.3). 

With Larmor gyrations ordered out of these equations, f s is independent of gyrophase. 
As a result, the pressure tensor is diagonal in a coordinate system defined by the parallel 
and perpendicular directions with respect to the magnetic field: 


P s = P±S 

(/ - 66) +p\\ s bb, 

(2.4) 

where / is the unit dyadic, and 



P\\s = n s T lls = 

J d 3 vm s (i>|| - U|| s ) 2 f s , 

(2.5) 

P-Ls = n s T ±s = 

f ]3 n 

/ d vm s -±fs 

(2.6) 


are the parallel and perpendicular pressures, respectively, of species s. By expanding 
Equation (2.3) in powers of (me/rrii) 1 / 2 and using the quasineutrality constraint, 


5> s n s =0, (2.7) 


it is straightforward to show that the parallel component of the electric field satisfies 



(V.P e \ 

1 r 

V en e y 

en e 


b-Vp\\ e 


(p. Le - P||e) V • b 


( 2 . 8 ) 


to leading order. f Equation (2.8) expresses parallel pressure balance for the (effectively 
massless) electron fluid, taking into account the forcing-out of large-pitch-angle particles 
(those with a = cos^ 1 (um/u) ~ 7t/ 2) from regions of increased magnetic-field strength. 
Indeed, in the absence of collisions, Equation (2.3) guarantees that each particle’s adi¬ 
abatic invariant p s = m s w^_/2B is identically preserved (i.e. f s remains constant along 
the phase-space trajectory Dp, s /Df = 0). The second term in parentheses on the left- 
hand side of Equation (2.3) is then straightforwardly interpreted as the mirror force, 
(m s Wj_/ 2)V •b = — fj.gb-'VB , which, by Equation (2.8), is compensated in the electron 
fluid by the parallel electric force, — eEy, and by the divergence of the parallel electron 
pressure, —(l/n e )V • (p|| e b). 


f Equation (2.8) may also be obtained directly from the electron force equation (Eq. 2.2 with 
s = e) after neglecting both the inertial terms on its left-hand side and the friction force on its 
right-hand side. 
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2.2. Particle Drifts and Centre-of-Mass Variables 

The same ordering as one uses to obtain a gyrophase-independent f s also guarantees that 
all species drift perpendicularly to the magnetic field with identical velocities. Indeed, 
the lowest-order contribution to the perpendicular electric field is E j_ = — ( u s /c ) X B , so 
that ttj_ s = u± = cE X B/B 2 . With this borne in mind, we interpret the D/D t operator 
in Equation (2.3) as measuring the rate of change of a quantity in a Lagrangian frame 
that is transported parallel to the magnetic field at velocity un and drifts perpendicular 
to the magnetic field at the E X B velocity. 

It then follows that the mean drift of any species relative to the centre-of-mass velocity 
n = E,. m s n s u s / Y, „ m s n s must be in the parallel direction, viz., u s = u + u'^ s b with 

u '\\s = y- s J d3v ( v \\- u \\)fs- (2-9) 

In centre-of-mass variables with p = m s n s , Equations (2.1) and (2.2) become 

+w v) p=-pV-u, (2.10) 



( 2 . 11 ) 


where p± = J2 S P± S and P\\ = are the total perpendicular and parallel pressures, 

respectively. The parallel drifts contribute to the total parallel pressure in a straightfor¬ 
ward way. The magnetic field satisfies the ideal induction equation 

f) R 

— = -cV xE = Vx(uxB), (2-12) 


i.e., the magnetic flux is frozen into a frame moving perpendicular to the magnetic field 
at the velocity u±. 


2.3. Reduced Ordering and Dimensionless Parameters 
We proceed by separating all fields into equilibrium values plus fluctuations: f s = fg s + 
6f s , B = B 0 z + SB , n s = n 0s + 6n s , p± s = p± 0s + <5p_i_ s , py s = p|| 0s + 5p|| s , and 
Up = wj| 0s + i5Mj| s . The mean magnetic field B 0 z and the equilibrium distribution function 
fo s are both taken to be spatially uniform (the slab limit). The asymptotic ordering 


*!L 

k± 

where 


SL 

f()s 


U± 

VA 


SB ± 

Bo 


at 

VA 


8B 1 

B 0 

VA = 


8n s 

Bo 


$P±s 

P±0s 


^Pj|£ 

P\\0s 


Su'n 


108 


(2.13) 


(2.14) 


y/Anpo 

is the Alfven speed, is then applied to “reduce” the equations so that they describe 
the evolution of anisotropic (fen <C k±) fluctuations whose parallel Alfven timescale and 
perpendicular non-linear timescale are of the same order, fcii^A ~ k±u± —the so-called 
critical-balance conjecture (Goldreich & Sridhar 1995), used here as an ordering assump¬ 
tion (S09). Such spatial anisotropy is both measured directly in the solar wind (e.g., 
Bieber et al. 1996; Horbury et al. 2008; Podesta 2009; Wicks et al. 2010; Chen et al. 
2011) and observed in numerical simulations of Alfvenic turbulence (e.g., Shebalin et al. 
1983; Oughton et al. 1994; Cho & Vishniac 2000; Maron & Goldreich 2001). The perpen¬ 
dicular perturbations are taken to be Alfvenic ( SB±/B 0 ~ u±/va) and the compressive 
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perturbations (6u\\, SB y, 5n a , Sp j_ s , <5py s , <5tiy s ) are ordered comparable to the Alfvenic 
ones, with the parallel and perpendicular thermal speeds of species s, 


12T\ 


^th||s 


||Q S 


and 


^th_l_s — 


2Tlo s 


(2.15) 


respectively, ordered comparable to the Alfven speed (i.e., /3y s , (3± s defined in Eq. 2.21 
ordered unity). We further assume that the characteristic frequency of the fluctuations 
to ~ knv\ ~ Vth|| i/L. This means that fast magnetosonic modes, for which ui ~ k±VA, 
are ordered out of our equations. Fast-wave fluctuations are rarely seen in the solar 
wind (Howes et al. 2012); observations of turbulence in the solar wind confirm that it is 
primarily Alfvenic (e.g., Belcher & Davis 1971) and that its compressive component is 
substantially pressure-balanced (Roberts 1990; Burlaga et al. 1990; Marsch & Tu 1993; 
Bavassano et al. 2004; see Eq. 2.26). We expect the same to hold true in the intracluster 
medium, where observationally inferred turbulent velocities are convincingly subsonic 
(e.g., Sanders & Fabian 2013; Zhuravleva et al. 2014). Indeed, our reduced ordering is 
consistent with a small sonic Mach number: Ma = u±/vthi ~ e(wA/^thi) 1. 

The density fluctuations of the various species are related to one another via quasineu¬ 
trality (Eq. 2.7): 


n 0s “ n 0s J 




(2.16) 


where c s == Z s no s /no e is the charge-weighted ratio of number densities; note that c e = — 1 
and c s = 0- The perturbed pressures are calculated via 


Sp\\s = J d^vi 


M w ll 


110s 


‘Sf. 


and 


Sp± s = / d 3 u )n s -T Sf s . 


(2.17) 


Perturbed parallel drifts may be obtained directly from taking first moments of the 
perturbed distribution functions, 


= — [ d 3 i> (u|| - uj, 0s ) 6f s - My with 
nos J 


Y* s m s / d 3 vv\\Sf s 
J2 S m sno s 


(2.18) 


rather than indirectly via the momentum equations (2.2) and (2.11). In other words, they 
are not independent quantities. 

The resulting set of equations has a number of dimensionless free parameters: 


71 


II Os 


71 


and 


11 Oe 


7j_os 

t -Ls ~ Tf, 

J J_0e 


(2.19) 


are the ratios of the parallel and perpendicular temperatures of species s to their respec¬ 
tive electron temperatures (note that 7]| e = r± e = 1); 


A s = ^ _ 1 

P||0s 


is the dimensionless pressure anisotropy of species s; and 

and 


. ^TPyoa 
P lb - 


B 2 0 


a ■ 87rp_L0 s 

P_Ls — 


B 2 0 


( 2 . 20 ) 


( 2 . 21 ) 


are the ratios of the parallel and perpendicular pressures of species s to the magnetic 
pressure. We use the shorthand /3y = ^) s /3y s an d P± = *• All of these quan¬ 

tities, evaluated in the equilibrium state, are taken to be order unity in the e expan¬ 
sion (subsidiary limits in, e.g., high and low /3y can be taken after the e expansion is 
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done; see §4.4). Likewise, we order uj| 0s ~ u t h||j for all species s. This precludes equi¬ 
librium parallel drifts in the electron fluid from entering into our equations, because 
uy ~ u t h|| e ~ •\/mi/TO e n t i 1 ||i i>th||i f° r the electrons (i.e., the random thermal motions 
of electrons are characterised by speeds much in excess of the ion thermal speed and, 
therefore, any parallel drifts in the equilibrium state). 

We will also make frequent use of the following compact notation: 


fll = _„ 2 

J Os ^th||s 


df 0s 


d(y\\ - u' 


II Os 


and 


f -L = _ v 2 dfos. 

hs th_U g w 2 ^ - 


( 2 . 22 ) 


which are dimensionless derivatives of a species’ equilibrium distribution with respect to 
the square of the parallel velocity (peculiar to the species equilibrium drift velocity) and 
perpendicular velocity (peculiar to the E X B drift velocity), respectively; and 


®fos 


P- LOs 
P||0s 


I Os 


fos, 


(2.23) 


which measures the velocity-space anisotropy of that distribution. These definitions are, 
of course, only useful insofar as particle collisions and scatterings are unable to maintain a 
Maxwellian distribution in the equilibrium state, for which /q s = = fo s and T)fo s = 0. 

To allow for this to be the case, we order the collision frequency va <C ew. This means 
that collisional relaxation of /o s towards a Maxwellian distribution occurs at higher orders 
than will be treated in this Paper, as does the heating of the background plasma due to 
collisional smoothing of the (secularly increasing) fine-scale structure in velocity space. 
For our purposes, the background equilibrium is thus stationary in time. 

While most applications to non-Maxwellian space and astrophysical plasmas make use 
of the bi-Maxwellian distribution function, 


fbi-M,s{v\\,W±) 



foil - «f,Q.) 2 

^thlls 


nv th±s 


exp 



w l \ 

V th±s ) 


(2.24) 


to describe the equilibrium distribution function of the plasma (for which /J s = f^~ s = /o s 
and D/os = A s f 0s ), we keep our derivation of KRMHD general with respect to the 
form of /os(v||,tujJ.f That being said, because of its widespread use by the space and 
astrophysical communities, we will often refer back to Equation (2.24) to present useful 
particular cases of our more general results. To ease application of our theory to the 
solar wind, in Appendix D we also specialise our equations for the case of a bi-kappa 
equilibrium distribution function. 

Finally, we caution that there are some situations in weakly collisional magnetised 
astrophysical plasmas for which the KRMHD ordering (2.13) is inapplicable. One such 
situation is when the equilibrium pressure anisotropy takes on values beyond the firehose 
stability boundary (Eq. 3.2), since k\\ is not much smaller than k± for the fastest-growing 
firehose modes. Likewise, beyond the mirror stability boundary (Eq. B13), the nonlin¬ 
ear evolution of the mirror instability involves the trapping of particles, a feature not 
accounted for in the KRMHD ordering (note however that the expected wavevector 


f Some restrictions on fo a are necessary in order to make sense of the generalised free-energy 
invariant of KRMHD, derived and discussed in Sections 4.1 and 5.1. Namely, /q s must be 
strictly positive throughout all of the parallel velocity space; i.e., fo s must decay monotonically 
away from uy = u(| 0s and not be too flat around that point. This restriction also eliminates the 
possibility that /o s is unstable to high-frequency bump-on-tail instabilities, which are outside the 
KRMHD ordering. This class of distribution functions covers all plausible distribution functions 
for the solar wind (Maxwellian, bi-Maxwellian, kappa, bi-kappa, etc.). 
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anisotropy k\\/k± <C 1 of these modes is captured). That being said, the majority of the 
solar wind does lie within these stability boundaries (e.g. Bale et al. 2009), and so our 
theory should be appropriate for describing Alfvenic turbulence in that part of parameter 
space. A second limitation of the KRMHD ordering is that there are times and places 
in both the solar wind and the magnetosheath where SB/Bo is not small (e.g. Alexan¬ 
drova et al. 2008), usually at the outer scale. However, as the free energy stored in the 
turbulent fluctuations cascades to smaller scales (in the inertial range and beyond), the 
fluctuations become smaller and more anisotropic. If anything, these fluctuations tend 
towards the KRMHD (or, more generally, gyrokinetic—see Appendix C) limit. 


2.4. Alfvenic Fluctuations 

To zeroth order in e, Equation (2.10) becomes Vx • u±_ = 0. Likewise, the divergence-free 
constraint on the magnetic field, V • B = 0, becomes Vx • SB± = 0. These simplifications 
allow the Alfvenic fluctuations to be expressed in terms of scalar stream (flux) functions: 

u± = z X Vl<1> and ^t=f= = z X VxT. (2.25) 

V^Trpo 

We substitute these expressions into Equations (2.11) and (2.12) and examine the result 
order-by-order in e. 

To lowest order, Equations (2.11) and (2.12) become, respectively, 


V± 



b^sbA 

4tt B q J 


= 0 , 


where the Poisson bracket 




(2.26) 

(2.27) 


{$,$} = i-(Vi$x Vit). (2.28) 

The first of these equations expresses perpendicular force balance, a result which will 
aid our description of the compressive fluctuations in Section 2.5. The second equation 
is identical to the induction equation in standard RMHD. At the next order, we derive 
an evolution equation for the stream function. Its simplest form is obtained by taking 
the ^-component of the curl of the force equation (2.11)— the vorticity equation—which 
gives 



For an isotropic equilibrium pressure (pxos = P||o SJ = 0) and no equilibrium parallel 
drifts (wj| 0s = 0), this reduces to the standard RMHD momentum equation. 

The nonlinearities in Equations (2.27) and (2.29) involving the magnetic field imply 
that Alfvenic fluctuations propagate along the locally deformed magnetic field rather 
than the uniform equilibrium field, and so the parallel and perpendicular directions do 
not strictly lie along the Cartesian axes defined by the guide field. Indeed, by introducing 
the Lagrangian operators 


d 

df 


d> 
dt 


+ u± • Vi 


d_ 

Ft 




(2.30) 
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(2.31) 


Equations (2.27) and (2.29) may be written compactly as 


d_ 
d t 



'A 


(2.32) 


(2.33) 


where we have introduced the effective Alfven speed 



1/2 


I’A* = 


(2.34) 


In systems where p±o — P\\o — w s no s wj 2 0s < 0, the speed at which deformations in 
the magnetic field are propagated is effectively reduced by the excess parallel pressure, 
which undermines the restoring force exerted by the tension of the magnetic-field lines. 
When p_lo ~P\\o ~ Y^ s m s n 0 su'^ 0s = —Bq/Ait, the magnetic tension is exactly balanced by 
the anisotropy in the distribution function and the plasma does not respond to (perpen¬ 
dicular) magnetic perturbations. Parallel drifts in the equilibrium distribution functions 
of different species make this criterion easier to satisfy by supplementing the parallel 
thermal pressure. For values of p±o ~P\\o below this threshold, the effective Alfven speed 
becomes imaginary and the stream and flux functions acquire a n/2 relative phase shift. 
The plasma then becomes firehose unstable, an effect that we discuss in Section 3.1. 
For now we caution that, if the equilibrium pressure anisotropy and parallel drifts make 
va*/va as small as e, the reduced ordering (2.13) is broken and the KRMHD theory, as 
derived here, becomes inapplicable. 

Equations (2.32) and (2.33) form a closed set, and so the Alfven-wave inertial-range 
cascade is completely decoupled from all other (compressive) types of perturbations (the 
Alfvenic cascade is further discussed in Section 3). While this result is usually derived 
starting from the collisional MHD limit, we have shown that the same holds true even 
for a collisionless plasma (as in S09) with arbitrary gyrotropic equilibrium distribution 
function. 


2.5. Compressive Fluctuations 


2.5.1. Parallel Electric Field 

To obtain the equations describing the density (Sn e ) and magnetic-field-strength (<5-B||) 
fluctuations, we return to the drift-kinetic equation (2.3). Applying the reduced ordering 
(2.13) and neglecting collisions, we have for the electron species 



(2.35) 


Further expanding Equation (2.35) in the small parameter (m e /77ii) 1,/2 removes the d/df 
term and the equilibrium electron drift (itj| 0e ) to lowest order. Dividing the result by uy 
and employing the 'D notation (Eq. 2.23) to group terms, we find 



(2.36) 
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Integrating over the velocity space leads to an expression for the parallel electric field, 


E I, = — 


1 


r 


h-V 


T\ 


||0e 


Oe 


( Sn e 


A 


le 


< 55,1 

Bo 


where we use the following notation: for integer £, 


I! _i_ 
°£e — 


1 

^Oe 




2 1 


& V ^th_Li 
= 1 for a bi-Maxwellian 


I Oe 


(2.37) 


(2.38) 


are dimensionless coefficients related to perpendicular moments of the parallel-differentiated 
equilibrium electron distribution function, and 

A^e = cj'e ~~ - 1 (2.39) 

P ||0e 

is the dimensionless pressure anisotropy (cf. Eq. 2.20) of the electrons weighted by those 
coefficients. For isotropic electrons, the parallel electric held is entirely related to fluc¬ 
tuations in the electron (and therefore ion) density; the corresponding (first) term in 
Equation (2.37) ultimately leads to the Landau damping of ion acoustic waves. When 
the equilibrium electron pressure is anisotropic, fluctuations in magnetic-held strength 
also contribute to the parallel electric held; this second term enforces quasineutrality in 
the face of preferential exclusion of large-pitch-angle electrons from regions of enhanced 
held strength. 


2.5.2. Pressure Perturbations 

With knowledge of the parallel electric held (Eq. 2.37), we can rewrite the perturbed 
electron distribution asf 




r 


Oe 


™0e 


A ^ 1 - 

le B 0 J foe 


w. 


SBn 


'th±e 


Bn 


S/ 0e 


(2.40) 


and compute the perturbed parallel and perpendicular electron pressures by taking the 
appropriate second moments (cf. Eq. 2.17): 


$P\\e 

P\\0e 

Sp±e 

P. LOe 


1 (5n e 

= —n- — + A 


cl \ n ° 


le 
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Oe 


^Oe 


le 


6B \\ 

Bo 

s Jh 

Bo 


-A« 


-2A 


Bo ’ 
SB\ { 


le 


Bo 


(2.41) 


(2.42) 


Note that, if the equilibrium distribution function of electrons is isotropic, no electron 
pressure anisotropy can be generated by the fluctuations and the electron huid remains 
isothermal along magnetic-field lines: 


b ■ VST e = 0. 


The latter occurs physically by rapid electron conduction along held lines. Deviations 
from isothermality in non-Maxwellian plasmas arise when electrons conserve their adia¬ 
batic invariant /i e in the presence of field-line compressions and rarefactions. Indeed, by 

f Technically, the perturbed electron distribution function can only be determined up to an 
additive unknown function whose parallel gradient vanishes. We have set this homogeneous 
solution to zero, a simplification which may be justified by assuming stochastic held lines (S09). 
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taking the third velocity moments of Equation (2.40), we see that the parallel flows of 
parallel and perpendicular electron heat, respectively 

Q\\ e = J d 3 vm e vf\5f e and Q ±e = j d 3 v m e v\\ ^ 8f e , 

satisfy b • VQ|| e = b • VQ ± e = 0. For bi-Maxwellian electrons, this translates into 

b-VST lje = 0 and b ■ V (dT ±e + A e {p e )6B^ = 0, (2.43) 

where (p e ) = m e <x e /2B 0 is the lowest-order contribution to the mean adiabatic invari¬ 
ant of the electrons. Equation (2.43) states that, while the parallel temperature of the 
electrons remains constant along field lines, the perpendicular temperature cannot do so 
without violating p e conservation (cf. eqs 39-40 of Snyder et al. 1997). 


2.5.3. Reduced Drift-Kinetic Equation 

It is often computationally convenient as well as physically illuminating to replace the 
perturbed distribution function 6f s by the function 


g s = Sf s 


w ± SB \\ f -L 

<l, Bo /0s 


(2.44) 


which is the perturbed distribution function if f 0s is taken to be a function of the exact 
adiabatic invariant p s = mw 2 j_/2B (rather than of mw\/2Bo)\ to wit, 


fs ~ f0s{v \\, p s ) = fs~ fos , 


m s w 


S | 


m s w 


S J 


- fs - / 0 s (U||,W_l) + 


2 B 0 2B 0 

w± 5B\\ dfo. 


-I 


2 B 0 dw_ l 


= Sf s - 


u th±. Bo 


fos — 9s- 


Indeed, using Equation (2.40) in Equation (2.44), we find that 


9e = 


C, 


Oe 


Sn t 


5B I, 


-nr —+ A le -^ - 


n-Oe 


Bn 


w 


2 5B II 


U t 2 h||e 


/Oe 


(2.45) 


(2.46) 


does not contain any derivatives of the equilibrium distribution function with respect to 
wj_. We will see that the same holds true for gi, whose evolution equation we now derive. 

The evolution equation for the perturbed ion distribution function is obtained by 
applying the reduced ordering (Eq. 2.13) to Equation (2.3) with s = i. Using Equation 
(2.37) for the parallel electric field, the reduced kinetic equation for the ions may be 
written in a compact form analogous to equation (145) of S09: 


Tt +v,b-V 


) Si + (i’ll - »||0i) 


b • V 


<5U|| 

d'Hli V n 0e le B 0 


y 0e 


w. 


5B ii 


<||i Bo 


fL = o. 


(2.47) 

Note that this equation does not contain any derivatives of the equilibrium distribution 
function with respect to wj_. In terms of gi , Equations (2.16), (2.18), and (2.26) become 
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■ 5u 'b = Y.cr / d3v ( v \\ - u 'woi)9i 

,lj 0i J 


(2.49) 


le 


Oe 




So 


E Ti. i c i 

71.nz 


d 3 t 




noi 


9i, 


L th_L?' 


(2.50) 


which reduce to equations (146)-(148) of S09 for a single-ion-species Maxwellian plasma. 

Equations (2.47)-(2.50) evolve the ion distribution function g *; the “slow-wave quanti¬ 
ties” it||, (JuG, and ABy; and the density fluctuations Sn e . All nonlinearities are contained 

in the d/d t and b ■ V Lagrangian operators, which include the Alfvenic quantities 4> and 
4/; these are determined separately and independently by Equations (2.27) and (2.29). 
Nonlinear scattering/mixing of slow waves and the entropy mode by the Alfvenic pertur¬ 
bations takes the form of passive advection of the distribution function gi. In other words, 
even when the equilibrium distribution function is non-Maxwellian and there are parallel 
drifts between the various species, the compressive fluctuations are passively transported 
by the Alfvenic fluctuations, a result that we have, thus, generalised from MHD (Lithwick 
& Goldreich 2001) and Maxwellian KRMHD (S09). The passive cascades of compressive 
fluctuations, as well as their kinetic damping and susceptibility to mirror instability in a 
pressure-anisotropic plasma are further discussed in Section 4. 


2.6. Summary 

The reduced theory derived here evolves 4 + N- lon unknown functions: 4>, 4>, 8B\\, 5n e , 
and gi for each of the N\ on different ionic species. The stream and flux functions, 4> and 
4/ respectively, are related to the fluid quantities (perpendicular velocity and magnetic- 
field perturbations) via Equation (2.25). They satisfy a closed set of equations, Equations 
(2.27)-(2.29), which describe the decoupled cascade of Alfvenic fluctuations whose phase 
speed is modified by pressure anisotropy and inter-species parallel drifts. In the collisional 
limit, they revert to the standard equations of RMHD. The density and magnetic-field- 
strength fluctuations (the “compressive” fluctuations, or the slow waves and the entropy 
mode in the collisional limit) require a kinetic description in terms of the ion distribution 
function g ,, which is evolved by the kinetic Equation (2.47). This kinetic equation itself 
contains 5n e and 8Bu , which are, in turn, calculated by taking velocity-space integrals 
of gi via Equations (2.48) and (2.50). The nonlinear evolution of gi, SB y, and 5n e is due 
solely to passive advection of g t by the Alfvenic turbulence, which mixes Sn e and 5B\\ in 
the direction transverse to the magnetic field. In the Lagrangian frame associated with 
the Alfvenic fluctuations, the compressive fluctuations obey a one-dimensional linear 
equation, which may be solved independently of the Alfvenic turbulence. 

Here we summarise our new set of equations: 


~dt 


va b • V4>, 


(2.51a) 


^vi$ = ^ A b-v4 ± vi^, 

at 


d t 


+ v l] b-V) gi + (v {l -UnoijS-V 


1 Zi 

cl T b 


Sn e . SBn 

-- + A ie - 

n 0 e Bo 


(2.51 b) 


w 2 ± ABy 

w th||i B ° 


4 = o. 

(2.51c) 
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d 3 v 


^0 i 


^th_L' 


-2* 

CjTj.iZi/ jO_Li 
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2A 2e - -^A le 
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Oe 
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u th±i 
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cL T±i , 


— hr- 2A - + +2ci E+ /- 

rj_i \ n II rdl / \ z —' c,-ti,-Zj/ du 


c, 


c, 


(2.51 d) 


(2.51e) 


(2.52) 


These equations reduce to equations (155)-(159) of S09 when the equilibrium distribution 
function is Maxwellian and only one ionic species is present (for which va* = ua, u \\oi = 0, 

c| e = !) T ||s = A ^e = o, u t h||s = v t hj_ s , /3|| s = /3_l s , and fl s = fo s = fos). 

It should be noted that Equations (2.51) are ideal, in that they are ignorant of any 
physics capable of dissipating the large gradients in phase space that will inevitably be 
produced as the turbulence cascades to smaller and smaller scales. While some of this 
fine-scale structure is regularised by finite-Larmor-radius effects, which are included in 
the non-Maxwellian gyrokinetic theory derived in Appendix C, numerical implementation 
of Equations (2.51) requires the addition of finite collisionality and resistivity. 


3. Alfvenic Fluctuations in the Inertial Range 

Having constructed a theoretical framework for the evolution of anisotropic kinetic 
turbulence in collisionless magnetised astrophysical plasmas, we now investigate its im¬ 
plications for the behaviour of Alfvenic fluctuations in the inertial range (dynamical equa¬ 
tions for these fluctuations were derived in Section 2.4). In this Section we demonstrate 
that pressure anisotropy and parallel drifts do not interfere with the nonlinear mixing 
of counter-propagating Alfvenic fluctuations. In doing so, we derive the two Alfvenic 
invariants that are independently conserved and cascaded by these interactions. The ex¬ 
planation of the physical content of these invariants is aided by the linear theory of Alfven 
waves, which we present in the next Section. 


3.1. Linear Theory of Alfvenic Fluctuations: Alfven Waves and Firehose Instability 
The linear theory of Alfvenic fluctuations in KRMHD can be readily obtained by dropping 
the nonlinear terms in Equations (2.51a) and (2.516) and adopting the solutions 4), 
4/ ~ exp(—i uit + i k • r). The resulting dispersion relation is simply 


W = ±fc|| 1 >A*, 

with eigenvectors satisfying $ = T^{va*/va)- When 

P. L0 -p\\o - 

s 


(3.1) 


(3.2) 


the phase speed of the Alfven wave becomes imaginary and the firehose instability results 
(Rosenbluth 1956; Chandrasekhar et al. 1958; Parker 1958; Vedenov & Sagdeev 1958). 
Physically, negative pressure anisotropies and/or parallel drifts reduce the elasticity of 
the magnetic-field lines, undermining the supplied restoring force necessary to propagate 
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the wave. In KRMHD, the fastest growth occurs at arbitrarily small parallel scales, 
with no small-scale regularisation accessible within the long-wavelength approximation, 
kpi <C 1, in which Equation (2.2) is derived. To obtain the fastest growing mode, finite- 
Larmor-radius effects must be taken into account (cf. Schekochihin et al. 2010). Direct 
calculation from the hot-plasma dispersion relation yields k^pi ~ |/3_l// 3|| — 1 + 2//3|| | 1//2 
for the parallel (k± = 0) firehose (Kennel & Sagdeev 1967; Davidson & Volk 1968) and 
kpi ~ 1 for the oblique firehose with k± ^ 0 (Yoon et al. 1993; Hellinger & Matsumoto 
2000 ). 


3.2. Elsasser Fields and Alfven Ratio 

The effective Alfven speed can be used to cast Equations (2.32) and (2.33) in a symmetric 
form via the introduction of the generalised Elsasser potentials , 

( ± =$±—'H, (3.3) 

va 


and the corresponding Elsasser fields 

z ± =zx VlC* 


v A * 5Bj_ 

«_L ±- 7 . 

V A V 47 T /90 


(3.4) 


The latter are a straightforward generalisation of the standard Elsasser (1950) variables 
to non-Maxwellian equilibria, for which va* ^ va- Combining Equations (2.27) and 
(2.29), one can show that the Elsasser potentials satisfy 


( d d \ 

\dt TVA *dz ) 


2 A± 


Vic 


-^({c + ,Vic} + {C , vic + } T Vi {C + ,C}). (3.5) 


Thus, the standard result that nonlinear interactions (“scatterings”) of Alfvenic fluctua¬ 
tions occur only between counter-propagating fluctuations (Kraichnan 1965) holds true 
for general (gyrotropic) distribution functions. What is modified by the non-Maxwellian 
nature of the distribution function is the amount of (perpendicular) magnetic fluctuations 
that comprise each of the Elsasser potentials. As va* —> 0, the magnetic fluctuations fail 
to propagate and the distinction between C + and C _ is no longer meaningful. Indeed, the 
very idea of critical balance that underpins the RMHD ordering kuVA ~ k±u± is based 
upon a causality argument: fluctuations cannot be correlated over a distance larger than 
that over which an Alfven wave propagates in a nonlinear interaction time. Significantly 
reducing the signal speed, with va*/va ~ e or smaller, interferes with this argument and 
breaks the reduced ordering used in this Paper. This is what will happen if the firehose 
threshold is approached. On the unstable side of the threshold, the firehose fluctuations 
that emerge are not anisotropic in the same way that Alfvenic, or more generally gy- 
rokinetic, fluctuations are: in fact, they have k\\ ~ k± (Yoon et al. 1993; Hellinger & 
Matsumoto 2000; Kunz et al. 2014). This is why the considerations in this Paper do not 
describe the turbulence on the unstable side of the firehose threshold. 

The fact that the Alfven ratio 


r a 




U th||a j 


(3.6) 


depends on the anisotropy inherent to the distribution function (cf. Barnes 1979) becomes 
testable in the solar wind, where measurements find that the energy in magnetic-field 
fluctuations exceeds the energy in the velocity fluctuations, ta < 1 (e.g., Belcher & 
Davis 1971; Matthaeus & Goldstein 1982; Bruno et al. 1985; Roberts et al. 1987; Tu 
et al. 1989; Marsch & Tu 1990; Grappin et al. 1991; Goldstein et al. 1995; Bavassano 
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et al. 1998; Podesta et al. 2007; Salem et al. 2009; Perri & Balogh 2010; Chen et al. 
2011; Borovsky 2012). This result is often interpreted in terms of “residual energy”, 
a r = (rA ^ 1)/(^A + 1), the difference between the energy in velocity and magnetic-field 
fluctuations that is believed to be an inherent feature of the turbulence itself (Pouquet 
et al. 1976; see Chen et al. 2013 and Wicks et al. 2013 for brief reviews of the relevant 
literature and contemporary analyses). While the observed scale-dependent component 
of the residual energy is likely to be intrinsic—recent theory predicts a oc k^ 2 residual- 
energy spectrum for both balanced and moderately imbalanced strong turbulence (e.g., 
Muller & Grappin 2005; Boldyrev et al. 2011, 2012)— the constant component can (at 
least partially) be attributed to non-MHD corrections to the Alfven speed due to pressure 
anisotropies and parallel drifts. Indeed, this interpretation is supported by Equation 
(3.6), and one may thus construct a more appropriate Alfven ratio by weighting the flux 
function by Va*/va (Belcher & Davis 1971; Matthaeus & Goldstein 1982). This was the 
route followed by Chen et al. (2013), who found that such a “kinetic normalisation”, 
equivalent to using 


rA* = 


va 

i’a* 'k 


(3.7) 


instead of ta, yields fluctuations that are closer to equipartition, with a mean residual 
energy of <r r * = (rA* — l)/(rA* + 1) = —0.19 and a mean Alfven ratio of ta* = 0.71 
(rather than cr r = —0.43 and ta = 0.40 using Eq. 3.6). Using the appropriate kinetic 
normalisation is thus essential when measuring quantities like the residual energy in 
plasmas with anisotropic distribution functions, such as the solar wind. 


3.3. Alfven-Wave Invariants 

Another standard result—that interactions between “+” and ” waves occur without 
exchanging energy—can also be shown to hold true in general. Multiplying Equation 
(3.5) by /OoC ± and integrating the result over space, we find 


dty AW 

df 


= 0, 


(3.8) 


where 



(3.9) 


are the independently conserved (free) energies of the forward- and backward-propagating 
Alfvenic fluctuations, respectively. Their sum, 


Waw = W+ w + W *w 

= \ /dVp 0 (|VACl 2 + |VAC-| 2 ) 

= \J d 3 r p 0 (JVa $| 2 + ^|Vi’k| 2 ^) 




w th||« J 



(3.10) 


is, of course, also conserved. If v a* (equivalently, the expression in the square brackets 
in Equation 3.10) is positive—i.e., if the plasma is firehose-stable— Waw is a positive- 
definite quantity measuring the total kinetic and potential energy stored in the Alfvenic 
fluctuations. To interpret Equation (3.10) when is driven negative by the pressure 
anisotropy and parallel drifts, we separate the various terms in the conservation law for 
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Waw as follows: 


d 
d t 



f Pouj 
V 2 


SBj_ \ — ^os^ios 


d 3 r 


A s 


fo/' 2 

ZU \\0s 

V t 2 h|| s 


d SB 2 ± 
dtl3f' 


(3.11) 


We then see that the terms on the right-hand side of this equation constitute a source for 
the kinetic and magnetic fluctuations on the left-hand side. Indeed, recent work on the 
nonlinear evolution of the firehose instability (e.g., Schekochihin et al. 20086; Rosin et al. 
2011; Kunz et al. 2014) has shown that the rate of relaxation of the pressure anisotropy 
A s is related to d5B\/dt and so one can, heuristically, interpret the right-hand side of 
Equation (3.11) as a velocity-space source of free energy multiplied by the rate at which 
fluctuations act to remove that source of free energy. In this case, the Alfvenic invariant 
(3.10) is minimised by growing fluctuations. 

Another way to interpret Equation (3.10) for a non-Maxwellian plasma is as follows. 
How close the equilibrium pressure-anisotropic distribution is to the firehose threshold 
has the effect of weighting the (free) energy associated with perpendicular magnetic 
perturbations. As the threshold is approached, bending field lines becomes energetically 
less demanding (dynamically, the negative-pressure-anisotropy stress cancels the tension 
force; see Eq. 2.11). As the threshold is crossed, Waw is no longer positive-definite and 
so can be conserved even if perturbations grow—which is indeed what happens. 

In Section 5.1, we show that Waw is part of a generalised free-energy invariant con¬ 
served and cascaded to small scales in phase space by the plasma turbulence. 


4. Compressive Fluctuations in the Inertial Range 

We now turn our focus to the behaviour of compressive fluctuations in the inertial 
range (dynamical equations for these fluctuations were derived in Section 2.5). In this 
Section, we show that the compressive fluctuations possess their own invariant, which 
has a natural interpretation when the equilibrium distribution function is cast in terms 
of the particle kinetic energy and adiabatic invariant. For a bi-Maxwellian plasma with a 
single ionic species, the inertial-range cascade of compressive fluctuations can be further 
split into three independent kinetic cascades. We derive the linear theory of compressive 
fluctuations in a bi-Maxwellian plasma (the general linear theory is given in Appendix B) 
and use it to demonstrate how pressure anisotropy affects the efficacy of linear parallel 
phase mixing and the partitioning of free energy amongst the various cascade channels. 


4.1. Compressive Invariant 

In Section 2.5.3, we derived the evolution and constraint equations for the perturbed ion 
distribution function gi, whose moments describe compressive fluctuations in a pressure- 
anisotropic plasma. The nonlinear evolution of these fluctuations is due solely to passive 
advection of gt. by the Alfvenic turbulence, which mixes Sn e and 5B\\ in the direction 
transverse to the local magnetic field. (Mathematically, this is a statement that, in the 
Lagrangian frame associated with the Alfvenic component of the turbulence, Equation 
(2.51c) is linear.) During this mixing, the compressive fluctuations satisfy an important 
conservation law, which we now derive. 

Multiplying Equation (2.47) by (T^/Z i )(c i /n 0i )(g i /f^) and integrating over the phase 
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space, we find that 



(4.1) 


T\\i w\ 5B\\ 


"th|h 


Bn 


= 0. 


On the other hand, multiplying Equation (2.47) by the term in the square brackets in 
Equation (4.1), integrating the result over phase space, and performing integration by 
parts gives 
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Using this expression in Equation (4.1), summing over ion species, and using Equations 
(2.48) and (2.50) to eliminate the resulting velocity-space integrals of g t produces the 
following conservation law: 



In the absence of interspecies drifts, Wcompr is an invariant conserved by Equations 
(2.47)-(2.50). The simpler version of fU compr that is conserved for the pressure-isotropic 
case (eq. 201 of S09) is related to the perturbed entropy SS s of the system (cf. Krommes 
& Hu 1994; Sugama et al. 1996; Howes et al. 2006; Schekochihin et al. 2008a). With 
interspecies parallel drifts, the right-hand side of Equation (4.3) constitutes a source or 
sink for this quantity. It is the work done by the fluctuating parallel electric held (Equa¬ 
tion 2.37) and by magnetic-mirror forces acting on the interspecies drifts, and represents 
the exchange of free energy between these drifts and the compressive fluctuations.) 

Adding and subtracting the phase-space integral of T^ 0e gl/^foe t° the right-hand side 
of Equation (4.4), with g e given by Equation (2.46), the compressive invariant may be 
re-written in the following compact form: 


Wcompr = / d d r 
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' d 3„W_(i +/ , x) Ei 

o7T 


2 /, 


Os 


(4.5) 


f See Appendix B.2 for a specific example of this physics, where we explicitly demonstrate 
that free energy can flow into or out of the interspecies drifts depending upon whether the 
system is unstable to an ion-acoustic (streaming) instability. 
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We now make a transformation analogous to that made in equation (149) of S09. Defining 


Sfs = g s + 


SB ii 


Wh|| s B o 


fL = Sfs + 


5B\\ 


^th-Ls 


Br 


■®/o« 


(4.6) 


and using Equations (2.46) and (2.50) to eliminate the resulting velocity-space integrals 
over g e and gt, respectively, we can rewrite Equation (4.5) in a particularly useful form: 
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(4.7) 


(4.8) 


is the extension of A 2e to arbitrary species s (cf. Eq. 2.39). This form of W CO mpr parallels 
the final expression in Equation (3.10) for the Alfvenic invariant: we have a quantity that 
has one interpretation if 1 — J2 S /3_LsA 2s > 0, namely that it is the generalised energy of 
the compressive fluctuations, and another if 1 — y~) /3 _l s A 2s < 0, in which case the final 
term in Equation (4.7) becomes a free-energy source for the mirror instability, for which 
the expression multiplying 5B^/8tt in Equation (4.7) is related to the stability parameter 
(see Section 4.4.2 and Appendix B.3); neglecting interspecies drifts, 


d_ 

dt 



T\\p s Sfs 
2 /<" 


Os 
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d SB\ 
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(4.9) 


In the unstable case, one can interpret dSB^/dt as the rate of relaxation of the pressure 
anisotropy as the mirror fluctuations grow (see, e.g., Schekochihin et al. 20086; Kunz 
et al. 2014; Riquelme et al. 2014; Rincon et al. 2015). In this case, Wcompr is minimised 
by growing fluctuations. 

Another interpretation of what happens when the stability threshold is approached and 
crossed is analogous to the one we offered at the end of Section 3.3 for a similar situation 
concerning the Alfvenic fluctuations. Within the mirror stability boundary, Wcompr is 
a positive-definite conserved free-energy-like quantity. As the system gets closer to the 
mirror threshold, it becomes energetically “cheaper” to produce magnetic compressions 
or rarefactions (<5Ry)—dynamically, this is due to the fact that the effect of positive 
Rlo — P||o is to reduce the magnetic pressure response (cf. Southwood & Kivelson 1993). 
Once the threshold is crossed, Wcompr is no longer positive definite and its conservation 
is compatible with the growth of SB n and Sf s (the mirror instability). 

The astute reader will recognise that the factor multiplying SB^ in the free-energy 
invariant, namely 1 — /3_l s A 2s , is not the exact mirror stability parameter, Equation 

(B 13). While the latter reduces to the former in the case of very high /3 or of cold electrons 
(for which the right-hand side of Eq. B 13 vanishes), in general there is a stabilising 
term due to the interaction of linearly resonant particles with the parallel electric field 
(i.e. Landau damping). This physics is contained inside the first term in the compressive 
invariant, proportional to d/ 2 . In order to see this, and to make better sense of the 
structure of Wcompr, we must understand the physical meaning of Sf s . 


4.2. Meaning of 5f s : (ny ,w±) vs. (e s , g s ) Coordinates 

Our decision to write fo s as a function of uy and w_ j_, while analytically convenient, is 
not the most natural choice for interpreting the compressive invariant Wcompr- Instead, 
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let us introduce the kinetic energy and adiabatic invariant of a particle of species s, given 
respectively by 



m s w 2 , 
^ = ^B- 


(4.10) 

(4.11) 


and rewrite the drift-kinetic equation (2.3) using e s and p s as our velocity-space coordi¬ 
nates (e.g., Hazeltine 1973): 


DA 

D t 


m s(n - uf| 0a ) 




d B , dB 

d i + U W° s ^ 


df s 

de s 



(4-12) 

where uy — uj| 0s = ±y/2 m s (e s — p s B). This is perhaps the most transparent form of 
KMHD: particles parallel-stream, E X iJ-drift, conserve p s , and change their kinetic 
energy by interacting with a parallel electric field and/or a changing magnetic field 
(mirroring) in the frame of the equilibrium species drift; the term — (Dw^/Dt) • b = 
u± • (Db/Df) is an inertial term having to do with the fact that the direction of the 
magnetic-field line changes as the particle streams along it and so the plane of the E X B 
drift tilts. 

When compared to the formulation of KMHD in (vinioj.) coordinates (Eq. 2.3), this 
formulation makes clear that there is a more general class of equilibrium solutions in 
a collisionless drift-kinetic plasma than those satisfying f s = fo s (v\\,w±), namely, f s = 
/os( £ siA t s)- Thus, what is referred to as the “equilibrium state” may, in fact, contain 
inhomogeneous fluctuations (e.g., Alfven waves), so long as the energy of each particle is 
conserved. As a result, there are pieces of 5f s in the (um,w;j_) formulation that may be 
absorbed into the equilibrium distribution function when fo s is taken to be a function e s 
and /j, s (e.g., fluctuations in magnetic-field strength that do not violate /^-conservation). 
To see that this is the case, we relate the two formulations via 
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SB I, 


^th_Ls Bo 


S/os. 


(4.13) 


Comparing Equations (4.6) and (4.13), it then becomes clear that the perturbed distri¬ 
bution function appearing in the compressive invariant Wc 0mpr (see Eq. 4.7) satisfies 

— ixA. 5B\\ ~ ~ 

Sfs = Sfs + —- D/os = Sfs + fos - fos = fs- fos, (4.14) 

v th±s Bo 


i.e., it is the perturbed distribution function if fo s is taken to be a function of (e s ,/Lt s ) 
instead of (u||, wj_). 

The meaning of the first entropy-like term in the compressive invariant Wcompr (see 
Eq. 4.7) is then readily apparent: it is the non-Alfvenic piece of the distribution function 
that represents changes in the kinetic energy of the particles due to interactions with the 
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compressive fluctuations. In it are contributions from Landau-resonant particles, whose 
energy is changed by the parallel electric and magnetic-mirror forces in such a way as to 
facilitate Landau (Appendix B.2) and Barnes (Appendix B.3) damping of ion-acoustic 
waves and slow modes. As long as the plasma stays within the mirror and streaming 
stability boundaries, the compressive invariant is positive-definite. 

In the next four sections (§§4.3-4.6), we show that the conservation of Wcompr repre¬ 
sents a turbulent cascade of compressive fluctuations from large to small scales in phase 
space. 


4.3. Parallel Kinetics: Two Decoupled Collisionless Cascades 

We begin by noting that, under our collisionless ordering, the u;_l dependence in Equation 
(2.51c) can be integrated out. We introduce two auxiliary functions, G n (y m) and G B (v m), 
defined so that 

6n e f SB\\ f 

-= / duy G n and —- = / duy G B , (4.15) 

noe J £> 0 J 

i.e., they contain all the w± integrals and species summations in the right linear com¬ 
bination as per Equations (2.51d) and (2.51e). Then Equation (2.51c) reduces to the 
following two coupled one-dimensional kinetic equations (cf. eqs 179 and 180 of S09): 


d Gn 

d t 


A v\\b • V G n 


E( v ii-“W)4(«ii)6-v 


^+ v \\b-VG b = || - «f| 0i )4(«||) b ■ V 
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\ nn (n i \ n B( n , \ 

X i ( u ll)— + A i hl)-5- 
noe i>0 , 



+ K“(v\\) 


(4.16a) 

SB \\ 

Bq \ 


(4.166) 


where the vm- and ion-species-dependent Xi coefficients, given in Appendix A, depend 
upon various perpendicular moments of the parallel-differentiated equilibrium ion distri¬ 
bution function: 

F « {vt) =w>L d —4(^) <4 ' i7) 

for integer L This coupled system of integral equations, compactly expressed as 
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has a simple physical interpretation. Kinetic fluctuations in each of the i’ = 1... A; on ionic 
species collectively excite compressive fluctuations in the magnetic-field strength (via 
perpendicular pressure balance), the electron density (via quasineutrality), and thereby 
the electric field (via parallel pressure balance in the massless electron fluid; see Eq. 2.37). 
These electromagnetic fields in turn feed back upon the kinetic fluctuations exhibited by 
each individual ionic species, in a way that is dictated by the various A j coefficients. This 
system of equations describing the evolution of the compressive fluctuations and their 
interactions with the particles is closed, signaling the fact that the compressive dynamics 
proceeds independently from that of the Alfvenic fluctuations. 

Because each ion species responds to the compressive fluctuations in a different way 
(Aj A i’ for i ^ i'), Equation (4.18) cannot be diagonalised in general. Moreover, the 
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mode structure excited in the parallel-velocity space by these fluctuations will generally 
be different than the mode structure present in the fluctuations themselves [in general, 
A(uy) ^ A(wm) for 'ey ^ i/,], and so this system can only be diagonalized if the \ co¬ 
efficients are independent of Dy. Both of these requirements for diagonalizing Equation 
(4.18) are guaranteed only for a plasma consisting of a single species of bi-Maxwellian 
ions, for which F|(uy) = F M (uy) = (1/v^thHi) exp(-u|/z; t 2 h y i ) for all integer L Then, 
assuming a bi-Maxwellian plasma and proceeding with the diagonalisation, we find that 
Equation (4.18) is equivalent to 


dG± .nS-VG^a^fe-V 


d£ 


A ± 


<Ml <?*(«) 


where 


G+ =G B + ^~ (l + —)c n and G~ = G n + — G B 

P_Li 


P-LOi 

are the eigenvectors; 
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T IK _|_ r llA ( P\\ 0 i Si 

TJft Zi J \P_L0i P_Li 


Zi P-LOi P_L 

are the (inverses of) the corresponding eigenvalues; and we have defined 


_ T \\i , r IK , P||o< 1 

=- h ^ H - ~a~ 

T_L i ^i PA-Oi PA_i 


2 / v 2 

T~\\i T\\i \ l P||Qi Q 


r±i 


\P-L0i /^_L? 


(4.19) 


(4.20) 


(4.21) 


(4.22) 


ft = 1 - P±A e , (4.23) 

Wi = l + ^A e (1 - A e ) , (4.24) 

P-LOz V 2 J 

the latter two expressions equating to unity for isotropic electrons. 

Equation (4.19), which reduces to equation (181) of S09 for a Maxwellian plasma, 
describes two decoupled kinetic cascades. These are the subject of Section 4.5, in which 
we derive the two collisionless invariants associated with G + and G~. But first, we 
specialise G ± and hA 1 for application to two different astrophysical systems representing 
two distinct parameter regimes and determine what their values imply for the evolution 
of compressive fluctuations in the linear regime. 


4.4. Linear Theory: Collisionless Damping and Mirror Instability 
To develop the linear theory for a two-component bi-Maxwellian plasma,) we Fourier 
transform Equation (4.19) in time (d/dt -» — iw) and space (b- V —>• ifcy), divide both 
sides by — i(w — feyt>y), and integrate over the parallel velocity. Dividing both sides of 
the resulting equation by / duyG 1 * 1 and using f dvy Fm(uii) = 1, we obtain the following 
dispersion relation: 

A± - 1 = f / dvy = SiZuiti), (4.25) 

k\\ J Uy - io/k y 

f The linear theory for a plasma with an arbitrary gyrotropic equilibrium distribution function 
and multiple ionic species is presented in Appendix B. 
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where = w/fci|Wth||» is the dimensionless phase speed and the (Maxwellian) plasma 
dispersion function (Fried & Conte 1961) 

1 e~ x2 

^m(0 = —j= / da;-(4.26) 

V 77 J -oo x - £ 

the integration being performed along the Landau contour. Formally, Equation (4.25) 
has an infinite number of solutions, most of which are strongly damped with damping 
rates Im(£j) ~ 1. A few of the more interesting solutions may be obtained analytically 
in the low- and high-beta limits. 


4.4.1. Low-Beta Limit: Solar Corona 

When /3 S <C 1 ~ A s , Equations (4.20) and (4.21) givef 

A--l~-^l + ^, G~ -G n + 7 ^ (l + Ai+Ae^J G b , (4.27) 

G+~G b . (4.28) 

P±0i P±i 

For the ” branch, we have Im(G) ~ 1, and so 


LO ~ 



(4.29) 


which is much smaller than the Alfvenic cascade rate k\\v a - For the “+” branch, consisting 
predominantly of fluctuations in magnetic-field strength, 


U! 


-i\k\\\v A \ P\\i 


In 


P\\0i 


P±0i /3j_i 


(4.30) 


up to logarithmically small corrections. This damping rate is slightly greater than that of 
the ” branch, though still much smaller than the Alfvenic cascade rate. Compressive 
fluctuations in a low-beta plasma are therefore weakly damped. Pressure anisotropies do 
not affect this conclusion. 


4.4.2. High-Beta Limit: Intracluster Medium 
When /3 S ~ 1/A S 1, we have 


A - — 1 ~ —2 [ + — 

, T±i Zii 


A+ - 1 ~ 


P||0i 


P±0i P±i 


(l-^x.A,), 


G~ ~ G„ 


G + ~ G b + \ — G n 

2 T\\i p±0i 


(4.31) 

(4.32) 


The ” branch corresponds to the density fluctuations and is strongly damped with 
Im(£j) ~ 1. In contrast, the damping rate of the “+” branch is small: it can be obtained 
by expanding Zm(&) = iy^ + which gives 


|fc||K'A P|| 0 i 

V^WiP±0i 


s 


(4.33) 


This expression generalises (in the limit ku/k± <C 1) for bi-Maxwellian distribution func¬ 
tions what is known in astrophysics as the Barnes (1966) damping and in plasma physics 


f In some regions of the solar atmosphere, this subsidiary expansion in low P s may conflict 
with the prior mass-ratio expansion if Pi ~ m e /mi. 
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as transit-time damping (Stix 1962). Particles that are almost at rest with respect to the 
slow wave (i.e. “Landau resonant” with <x> ~ feiiUM) are subject to the action of the mirror 
force associated with the magnetic compressions in the wave. Since, for a monotonically 
decreasing distribution function (/^ > 0), there are more particles with dm < w/k\\ than 
with U|| > w/fcy, the energy exchange between resonant particles and the wave leads to a 
net gain (loss) of energy by the particles (wave). Put differently, the only way to maintain 
perpendicular pressure balance for a slow wave is to increase the energy of the resonant 
particles (referred to as “betatron acceleration”, due to the last term in square brackets 
in Eq. 4.12) at the expense of the wave energy. The result is wave damping. 

Because G n is strongly damped, the fluctuations that are damped at the rate (4.33) 
are predominantly of the magnetic-field strength. For large ion beta, the damping rate 
is a small fraction ~l/dE of the Alfvenic cascade rate, becoming even smaller in a 
plasma exhibiting positive pressure anisotropy. This reduction is due to the proportional 
increase (for A s > 0) in the number of large-pitch-angle particles in the magnetic troughs 
(SB || < 0). This inflates the field lines (in order to maintain perpendicular pressure 
balance), thereby (partially) offsetting the damping of the field-strength fluctuations. If 
the concentration of these particles leads to more perpendicular pressure than can be 
stably balanced by the magnetic pressure, the troughs must grow deeper to compensate. 
This process runs away as the resonant particles in the deepening troughs lose energy 
at a rate pdB/dt ~ fi'ySB^ (“betatron deceleration”). This is the mirror instability 
(e.g. Southwood & Kivelson 1993). It is the inevitable outcome of trying to maintain 
perpendicular pressure balance in the midst of an effectively negative magnetic pressure 
(see Appendix B and, in particular, Eq. B6). 


4.5. Collisionless Invariants and Phase Mixing 

The linear theory elucidated, we now return to the principal result of Section 4.3: that, 
for a plasma with a single species of bi-Maxwellian ions, the compressive fluctuations 
can be decomposed into two decoupled kinetic cascades. Here we show that these two 
cascades independently obey their own conservation laws. 

If we multiply Equation (4.19) by G' ± (t!||)/i 7 M('C||), integrate over space and parallel 
velocity, and perform integration by parts on the right-hand side, we find that 

Ttj iSr J dv 11 2^J) = - A ^ / d 3 "- (/ d »l G±) 6 • v / d», «, G±. (4.34) 

On the other hand, integrating Equation (4.19) over parallel velocity gives 

^ J du,| G ± = —fo • V J d«,| v,| G ± , (4.35) 


because f duy t>i|i i M(' t, ||) = 0. Using this to replace the final term on the right-hand side 
of Equation (4.34), we find 


where the two invariants are 


dW ± 

r compr _ 

dt = ’ 


Ti-± _ n ojT\\oi 
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d 3 r 




(4.36) 


(4.37) 


For these invariants to be conserved, the first term in the square bracket must grow to 
compensate for the decay of the second due to collisionless damping. To see how this 
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arrangement proceeds, it is useful to split 


G ± = G ± + / du||G ± , 


noting that J dvn G ± = 0 since f dwn jFm(vi|) = 1- Then 


(4.38) 


± _ no,T l]0l 

^'compr c\ 


d 3 r 





(4.39) 


For a Maxwellian equilibrium distribution function, the two invariants W± mpl are guar¬ 
anteed to be positive definite because A + > 1 and A - < 0 (see Eq. 4.21). Collisionless 
damping leads to exponential decay of the density and magnetic-field-strength fluctua¬ 
tions, or, equivalently, of f dt>| G ± , while conserving W± mpT . This means that damping is 
a redistribution of the conserved quantity W,~ mpr : the first term grows to compensate for 
the decay of the second. This is a manifestation of linear parallel phase mixing (Landau 
1946; Hammett et al. 1992; Krommes & Hu 1994; Krommes 1999; Watanabe & Sugama 
2004): free energy passes from the low (density) moment of the distribution function 
to higher moments (contained in G ). As time goes on, the latter part of the solution 
becomes increasingly oscillatory in un (G^ oc e _lfe u*ll i , the so-called ballistic response), 
representing the development of finer structure in the parallel-velocity space (see §6.2.4 
of S09 for further discussion). 

What is novel for anisotropic distribution functions is what happens near the mirror 
threshold A + = 1 (cf. Eq. 4.32). When the perpendicular pressure is greater than the 
parallel pressure, the Barnes damping rate is reduced (Eq. 4.33), and so the free energy 
is transferred from the electromagnetic fluctuations to the ballistic kinetic fluctuations at 
a slower rate. As we explained in our discussion Equation (4.7), it is easier to compress 
magnetic-field lines (5B y ^ 0) when p lo > P\\o, so the energetic cost of perturbing the 
magnetic-field strength is relatively small; this is the physical origin of the small factor 
(1 — 1/A + ) multiplying the second term in Equation (4.39). In other words, as the mirror 
threshold is approached, the generation of fine-scale structure in velocity space by phase 
mixing is achieved by damping increasingly large magnetic-field-strength fluctuations 
at an ever decreasing rate. When A + < 1 and the plasma is driven mirror-unstable, 
the second term in Equation (4.39) becomes increasingly negative as the magnetic-field- 
strength fluctuations grow. Consequently, the first term in Equation (4.39) must grow 
increasingly positive to compensate, which corresponds to the production of fine-scale 
structure in velocity space. This structure is caused by the fraction of particles that 
are linearly resonant with (and nonlinearly trapped by) the unstable mirror mode (e.g. 
Southwood & Kivelson 1993). 


4.6. Wcompr Revisited: Compressive Phase-Space Cascade 
We now show that, for a bi-Maxwellian plasma with a single ion species, the compres¬ 
sive invariant Wcompr obtained in Section 4.1 incorporates the two invariants derived in 
Section 4.5. We begin by expressing the density and magnetic-field-strength fluctuations 
in terms of G ± : 

wci O'* /*’" Av 11G+ ) • <440) 
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where <r,; and Wi are defined in Equations (4.22) and (4.24) andf 
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To express g t in terms of G ± , we decompose gi as follows (cf. §6.2.5 of S09): 
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(4.42) 


(4.43) 


£=0 


where x = w^_/v^ h±i and the Laguerre polynomials Lt{x) = {e x / £\){<r / dx^) are x . Be¬ 
cause Laguerre polynomials are orthogonal, we have 
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(4.44) 


where the expansion coefficients are determined via the Laguerre transform: 

/»00 

G^(uy)= / dxe~ x L e (x)g(x,v\\). (4.45) 

Jo 

Since L 0 = 1 and L\ = 1 — x, Equations (2.48) and (2.50) can be used to write the £ = 0 
and £ = 1 expansion coefficients as linear combinations of G n and Gb- Using Equations 
(4.20) and (4.21) to replace G n and Gb by suitable combinations of A + G + and A ~G~, 
we find after some straightforward but tedious algebra that 
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(4.46) 


(4.47) 


Substituting Equation (4.43) into the ion kinetic equation (2.51c), we see that all higher- 
order expansion coefficients satisfy a simple homogeneous equation: 


d Gt 
df 


> 1. 


Thus, the distribution function can be explicitly written in terms of G ± : 
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(4.48) 


(4.49) 


where g t comprises all Gi with £ > 1. In other words, g t is a passively mixed, undamped, 
ballistic-type mode that contributes to neither density nor magnetic-field strength: 


d. 9i 
d t 


+ v\\b-Vgi = 0, d 3 vgi=0, 


d 3 u g t = 0; 
Uh_L i 


(4.50) 


it is the homogeneous solution of Equation (2.51c). 

Equipped with expressions for Go (Eq. 4.46), G\ (Eq. 4.47), Sn e /no e (Eq. 4.40), and 

f There is a factor of 2 missing in the definition of k in S09 (their eq. 204). This error affects 
their equation (213). 





















26 M. W. Kunz, A. A. Schekochihin, C. H. K. Chen, I. G. Abel, & S. C. Cowley 


5B\\/Bq (Eq. 4.41) written in terms of G , we now substitute these into our expression 
for the compressive invariant (Eq. 4.5 with Eq. 4.44) to find that 


Wcompr = d 3 r 
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Thus, for a bi-Maxwellian plasma, the generalised invariant for compressive fluctuations 
splits into three independently cascading parts: W± mpr associated with the density and 
magnetic-field-strength fluctuations and a purely kinetic part given by the first term in 
Equation (4.51): 

W Th = J d 3 r j d 3 v T ^B. (4.52) 

All three cascade channels lead to small perpendicular spatial scales via passive mixing 
by the Alfvenic turbulence and to small scales in V|| via the linear parallel phase mixing, 
the rates of mixing being functions of the velocity-space anisotropy of the equilibrium 
distribution function. 


5. Conclusions 


5.1. Generalised Free-Energy Cascade in a Pressure-Anisotropic Plasma 

Assembling the results of Sections 3 and 4, we now arrive at the central unifying concept 
of this Paper. The Alfvenic invariants (Eq. 3.10) and the compressive invariant (Eq. 4.7) 
together make up the generalised free energy , 


W = TT+ W + W^w + Wcoznpr 
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(5.1) 


which is the quantity conserved by Equations (2.51) in the absence of equilibrium in¬ 
terspecies drifts and cascaded to small scales in phase space across the inertial range of 
KRMHD turbulence (analogous to the energy cascade in fluid or MHD turbulence). It 
contains (in the order of appearance in Eq. 5.1) the perturbed entropy of the system in 
the frame of the Alfvenic fluctuations, the energy associated with the E X B motion, the 
energy carried by the magnetic fluctuations, and terms arising from the exchange of free 
energy between the magnetic fluctuations and the equilibrium pressure anisotropy and 
parallel interspecies drifts. Just as shown by S09 for a two-species Maxwellian plasma, the 
inertial-range kinetic cascade can generally be split into three independent cascades of 
the generalised Alfvenic and compressive-fluctuation energies: W^w, and W com pr- 

In Section 4.6, we showed that, for a single-ion-species bi-Maxwellian plasma, Wc 0mpr can 
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be further decomposed into three independently cascading parts: W ,+ mpr , W c “ m pr> an d 
Wg r In all of these cascades, what is affected by pressure anisotropies and interspecies 
drifts are the amount of magnetic-field fluctuations associated with each of the invari¬ 
ants (§§3.3, 4.3) and the rate at which linear parallel phase mixing generates small-scale 
structure in velocity space (§4.4). 

Because we have ordered collisions out of our equations, the collisionless invariant given 
by Equation (5.1) is just one of an infinite number of invariants of the system. We place 
special emphasis on this particular invariant, W, for two reasons. 

First, it neatly reduces to the generalised invariant for a (collisional) Maxwellian plasma 
(cf. eqs. 74 and 153 of S09): 


W ->• 



Pouj |AB| 2 \ 

2 8tt J ’ 


which is variously referred to as the generalised grand canonical potential (Hallatschek 
2004) or free energy (Fowler 1968; Scott 2010) because of its similarity to the Helmholtz 
free energy A = — T 0s SS s + SU , where SU is the potential energy stored in the 
fluctuations and SS s is the entropy associated with the perturbed distribution function. 

Secondly, it encodes rather neatly in a thermodynamical context the main (linear and 
nonlinear) physical effect associated with the presence of firehose (§3.1) and mirror (§4.4) 
instabilities. The firehose and mirror (in)stability parameters appear in Equation (5.1) 
as prefactors of the perpendicular and parallel magnetic energies, respectively. As either 
of these thresholds is approached, these pre-factors get smaller and thus the energetic 
weight associated with the corresponding type of magnetic fluctuation gets smaller. At 
the firehose threshold, bending magnetic-field lines is free; at the mirror threshold, com¬ 
pressing them is free. Beyond these thresholds, it becomes energetically profitable to 
grow magnetic fluctuations. 

As long as the plasma stays within the firehose and mirror stability boundaries, the 
generalised free-energy invariant W is positive-definite and so describes a turbulent cas¬ 
cade from large to small scales (as well as into phase space). In this case, linear stability 
implies nonlinear stability. As these thresholds are crossed, the reduced ordering un¬ 
derpinning our equations breaks down, since there are no restoring forces to keep such 
fluctuations small. Thus, Equation (5.1) can only be meaningfully interpreted as a free- 
energy invariant when the plasma is stable. 

In the presence of equilibrium interspecies drifts, the conservation law for the gener¬ 
alised free energy acquires source/sink terms: 


d W 
d t 
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(5.2) 


which correspond to the change in the free energy due to the work done on the system 
by the fluctuating parallel electric (Eq. 2.37) and by magnetic-mirror forces acting on 
the equilibrium parallel drifts. 


5.2. Quantitative Details Matter 

Within the stability boundaries imposed by the firehose, mirror, and streaming instabil¬ 
ities, the theory presented here is analogous to that presented in S09 for a single-ion- 
species Maxwellian plasma, with many of the differences amounting to a somewhat mun¬ 
dane change of multiplicative coefficients. But this in itself is interesting, for it demon¬ 
strates that many of the salient qualitative features of the gyrokinetic and KRMHD 
theories of astrophysical turbulence are robust with respect to deviations from velocity- 
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space isotropy. On a more practical quantitative level, these otherwise benign coefficients 
do affect, at the order-unity level, the relationships between various fluctuating fields 
and partitioning of free energy that lie at the heart of predictive theories of solar-wind 
turbulence. In light of the ever-increasing scrutiny placed upon theories of Alfvenic tur¬ 
bulence by the wealth of data from the solar wind, as well as the astrophysical importance 
of knowing the proportion of turbulent energy that is distributed between the ion and 
electron populations, such details matter. 

The equations derived and discussed herein may be readily incorporated into existing 
(Maxwellian) KRMHD and gyrokinetic numerical codes. Such an advancement is now 
underway. 

Support for M. W. K. was provided by NSF through Max-Planck/Princeton Center for 
Plasma Physics Grant No. PHY-1144374 and by NASA through Einstein Postdoctoral 
Fellowship Grant No. PF1-120084, the latter issued by the Chandra X-ray Observatory 
Center, which is operated by the Smithsonian Astrophysical Observatory for and on 
behalf of NASA under Contract No. NAS8-03060. Support for C. H. K. Chen was pro¬ 
vided by an Imperial College Junior Research Fellowship. The authors thank Michael 
Barnes, Ron Davidson, Greg Hammett, Anjor Kanekar, John Krommes, Felix Parra, Ja¬ 
son TenBarge, and George Wilkie for useful conversations. The completion of this work 
was facilitated by the generous hospitality and material support provided by the Inter¬ 
national Space Science Institute (ISSI) in Bern to ISSI/ISSI-Beijing Team 304 and by 
Merton College, Oxford. 


Appendix A. Definitions of Aj-coefficients in Equation (4.16) 

In Section 4.3, we derived Equation (4.18) describing the parallel kinetics of the com¬ 
pressive fluctuations, which explicitly shows how kinetic fluctuations evolve under the 
influence of the density and magnetic-field-strength fluctuations they excite. This Equa¬ 
tion involves 4Aj on coupling coefficients, which we give here in terms of the usual plasma 
parameters and various perpendicular moments of the parallel-differentiated equilibrium 
ion distribution function, denoted (see Eq. 4.17): 
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where the bracket [... ] is given by Equation (2.52). 
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Appendix B. Linear KRMHD Theory for Arbitrary fo s 

In Section 4.4 we derived the linear theory of compressive fluctuations under the as¬ 
sumption that the plasma contained only a single species of bi-Maxwellian ions. This 
simplification was necessary in order to diagonalize the ion kinetic Equation (4.16) and 
thereby decompose the compressive invariant Wcompr into its three independently cascad¬ 
ing parts. In this Appendix, the linear theory for the compressive fluctuations is derived 
for a plasma containing multiple ionic species with arbitrary Jq s . 


B.l. Linear Dispersion Relation 

To obtain the linear dispersion relation governing the compressive fluctuations, we return 
to Equation (2.51c). Dropping the nonlinear terms and Fourier transforming in time 
(d/dt —► —iw) and space (b • V —> ifcy), we find 


9ik 


1 f ek . a 

+ le ^J 


w\ 6B\\ k U|| m ||Qj || 
v th\\i B ° \ u ll - w / fc ll ^ 


(Bl) 


Computing the required moments of this equation introduces poles where ui = k»vu, 
which correspond to wave-particle resonances and generically result in the growth or 
decay of compressive fluctuations. Anticipating this, we define the ion coefficients 
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which engender suitable generalisations of the plasma dispersion function for non-Maxwellian 
distributions: e.g., 


Cu = 1 + Ci^M(Ci) for a bi-Maxwellian, where & 


u ~ k W U \\0z 

l^iiKhii* 


(B 3) 


Note that these ion coefficients resemble the corresponding electron ones (cf. Eq. 2.38) 
in the limit ui, fc||itj| 0e -C fe||W|j, By analogy with Equation (2.39) for the electron fluid, we 
also define the dimensionless pressure anisotropy of the ions appropriately weighted by 

(j\\. 

A* = 4^-i- ( B4 ) 

P||Oi 

Note that the effect of equilibrium parallel drifts is implicit in the modified pressure 
anisotropy A&, as they enter through the dependence of the C H coefficients. 

With these definitions in hand, we proceed as follows. Taking the zeroth and miw \/2 
moments of Equation (B 1) gives, respectively, 
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Equation (B6), along with Equation (2.42) for the perturbed perpendicular electron 
pressure, allows the pressure-balance relation (Eq. 2.26) in the linear regime to be written 
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as 


Cia Z s (3± s {Sn ek . SB l{k \ 

V^o7 le ^J 



(B7) 


Combining this with the linearised quasineutrality equation (Eq. 2.16 with Eq. B5), 
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ClZs_ 



we obtain the KRMHD dispersion relation 



(B 8) 


(B9) 


The left-hand side of this equation exhibits the usual modification of the plasma beta pa¬ 
rameter by the pressure anisotropy of each species. For a single-ion-species bi-Maxwellian 
plasma, Equation (B9) reduces to Equation (4.25), which most notably contains the 
pressure-anisotropic version of the Barnes damping (Eq. 4.33). We now specialise Equa¬ 
tion (B 9) for two interesting cases. 


B.2. Landau Damping and Ion Acoustic Instability 
Let us first consider the /3 m —>• 0 electrostatic limit of Equation (B9): 
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We specialise Equation (B 10) for a Maxwellian plasma (cf. Eqs 2.38 and B 3) consisting 
of massless electrons, cold ions (|£j| 1 ), and a small population (c a -C 1 ) of drifting 

( u || 0 a < Uh||a) hot alpha particles (|f a | < 1). Expanding Z M (&) - -1/& - 1/2and 
■^M(Ca) — ii/r an< l assuming that the decay/growth rate 7 is much smaller than the 
real part of the frequency w r , we find 
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(B 11) 

where v s is the ion sound speed. For w r > fc||tin 0a , this is simply Landau damping of an ion- 
acoustic wave, occurring at a reduced rate due to the flattening of the total distribution 
function in the vicinity of the alpha-particle drift velocity. For w r < k^u'n Qa , Equation 
(B 11) represents a sort of ion-acoustic instability, driven by the differential streaming of 
alpha particles (rather than of electrons, which is the standard case; e.g. Davidson 1983). 
(Note that the assumption 7 < w r is satisfied because c a <C 1.) Outside of the limits 
taken, this branch of the dispersion relation is heavily damped (§4.4.1). 

In Section 4.1, we stated that the right-hand side of the evolution equation for the 
compressive invariant (Eq. 4.4) represents the work done by the fluctuating parallel 
electric field and by magnetic-mirror forces acting on the interspecies drifts, and thus 
is related to the rate at which free energy is exchanged between these drifts and the 
compressive fluctuations. Here we demonstrate this explicitly for the electrostatic case 
investigated in this section. Using Equation (B 5) to express the density perturbation 
of the alphas in terms of the density perturbation of the electrons, expanding Cq 0 ~ 
1Ti^/and using Equation (B 11) for the real and imaginary parts of £ a , we find that 
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our free-energy equation (Eq. 5.2) becomes 
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(B12) 


When the plasma is streaming-unstable, namely ui T < fc||wj| 0a j the right-hand side of this 
equation becomes positive, and free energy is extracted from the interspecies drifts and 
put into the compressive fluctuations. 


B.3. Barnes Damping and Mirror Instability 


Next, we treat the linear theory for Barnes damping and the mirror instability. Two 
results are sought: (i) the general mirror stability threshold for a single-ion-species non- 
Maxwellian plasma; and ( ii ) the decay/growth rate for a bi-Maxwellian plasma consisting 
of massless electrons, hot ions (|£j| <C 1), and a small (c a <C 1) population of drifting 
( u '\\ 0 a < t>th||a) hot alpha particles (|£ a | < 1). 

For («), we note that the transition from stability to instability proceeds through w —> 0 
(this is not generally true in the case with particle drifts, as we show below). In this limit, 
the coefficients 
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for both the ion and electron species, and the stability criterion may be read off directly 
from Equation (B 9) : 



(B 13) 


This is in agreement with existing expressions in the literature (e.g. Pokhotelov et al. 
2002; Hellinger 2007, eq. 34). 

For (ii), we expand Zm(^) — Zm(£o) — Iv/tt in Equation (B9) and note that £ a can 
be expressed in terms of 
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Solving for the real and imaginary parts of the frequency, we find 
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(B 156) 


Thus, the usual decay/growth rate (Eq. 4.33) is modified by the pressure anisotropy of 
the alpha particles and reduced by the effect of the hot alphas on the Landau resonance. 
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The mode acquires a real part proportional to the alpha-particle drift so that, even at 
marginal stability, w/0. 


Appendix C. Derivation of KRMHD from Pressure-Anisotropic 
Gyrokinetics 

Heretofore, we have worked under the assumption that both w/f2 s and kp s are asymp¬ 
totically small to all orders in e, where fl s = q s Bo/m s c is the gyrofrequency and p s = 
^th_Ls/Ps is the gyroradius of species s. In this Appendix, we relax this assumption and al¬ 
low for fluctuations with k±p s ~ 1 and to ~ ef2 S . The resulting set of nonlinear equations 
generalises Maxwellian slab gyrokinetics (Howes et al. 2006; S09) to arbitrary equilibrium 
distribution function /q s . Some of these results have been obtained before—notably, the 
nonlinear gyrokinetic equation for arbitrary /q s was derived already by Frieman & Chen 
(1982). Our restriction to slab geometry removes many of the complications introduced 
by those authors’ applications to axisymmetric tokamaks and makes the theory rather 
easier to grasp, which hopefully gives this Appendix some pedagogical value. We will also 
use these equations to show that the generalisation of KRMHD derived in this paper can 
be rigorously obtained from a gyrokinetic description by taking the k±p s <C 1 limit. 


C.l. Basic Equations and Gyrokinetic Ordering 
We begin with the kinetic Vlasov-Landau (or Boltzmann) equation 
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The electric and magnetic fields are expressed in terms of scalar and vector potentials: 

and B = B 0 z + V X A, (C 2) 
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where V- A = 0 (the Coulomb gauge). In the non-relativistic limit, these fields satisfy 
the plasma quasineutrality constraint (which follows from the Poisson equation to lowest 
order in k 2 Ajy where Ad is the Debye length), 


0 = ^ q s n s =$> / d 3 v f s , 


and the pre-Maxwell version of Ampere’s law, 


- V 2 A = — j = —^q s f d 3 vvf s . 


(C 3) 
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In this Appendix, and in contrast with the main text, we work with the velocity-space 
coordinate v ±, i.e., the full velocity perpendicular to the magnetic-field direction, rather 
than with the perpendicular velocity w± = v± — u± peculiar to the E X B flow 
Equations (C1)-(C4) are reduced by expanding 


fs — fos + Sfi a + 5f2s + ■ ■ •, 


(C5) 


where the subscript indicates the order in e, and by adopting the ordering of Section 2.3 
along with k±p s ~ 1 and u> ~ efl s : 


w p s 

h 

u± 

6B ± 

U\\ 

SB\\ 

n s l 

k_ l 

VA 

B 0 

VA 

B 0 

k±Ps 

~ d\\s 

~ d±s 

~A S 

~ T ||s ' 

~ T_Ls 


(C 6) 




33 


Kinetic Turbulence in Pressure-Anisotropic Plasmas 

We remind the reader that the collision frequency Va -C e 2 fl s , allowing for non-Maxwellian 
fos (cf. §A2.2 of Howes et al. 2006). Also, because the electron-ion mass-ratio expansion 
is not performed at the outset of the gyrokinetic derivation {m e /mi <C 1 is taken as a 
subsidiary expansion in §C.8), we must allow for the possibility of parallel electron drifts 
in the equilibrium state (wj| 0e ~ u t h||e)- 

The formal expansion of Equations (C1)-(C4) that results is worked out order by 
order. We begin with Equation (C 1), ordered relative to ujfo s - 


C.2. Gyrokinetic Equation 
C.2.1. Minus-First Order, 0( 1/e) 

The largest term in Equation (C 1) corresponds to Larmor motion of the equilibrium 
distribution about the uniform guide field: 

-d s z^rx|^j= 0 . (C 7 ) 

Decomposing the particle velocity in terms of the parallel velocity un, the perpendicular 
velocity v±, and the gyrophase angle i9, 

v = vnz + uxfcos^a: + sini?y), (C8) 

we find that Equation (C7) takes on the simple form 

- Qs lW = 0 - (C9) 

The equilibrium distribution function is thus independent of gyrophase (gyrotropic): 

fos = fos(v\\,v±,t). 


C.2.2. Zeroth Order, 0(1) 

Proceeding to next order and decomposing the velocity into its parallel and perpen¬ 
dicular parts (Eq. C8), Equation (Cl) becomes 
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This equation is simplified by employing the shorthand (cf. Eqns 2.22 and 2.23) 
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where the order in e of each term is indicated. The resulting lowest-order equation for 
Sfla, 
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admits a homogeneous solution and a particular solution. 
The homogeneous solution h s satisfies 
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where we have transformed the i9 derivative taken at constant position r to one taken at 
constant guiding centre—the centre of the ring orbit that the particle follows in a strong 
guide field: 

R s =r+ 1 ^. (C15) 

Thus, h s is independent of the gyrophase angle at constant guiding centre R s (but not 
at constant position r): 

h s =h s (t,R s ,v ||,ul). (C 16) 

It represents the response of charged rings to the perturbed fields, and is thus referred 
to as the gyrokinetic response. 

The particular solution of Equation (C 13)— the so-called adiabatic, or “Boltzmann”, 
response—is given by 
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It arises from the evolution of f 0s under the influence of the perturbed electromagnetic 
fields, a fact that is most clearly demonstrated by transforming Equation (C 17) to (£ a ,p s ) 
coordinates (see §C.4). It is instructive to note that the combination 


V's = V ~ u '\\Qs A \\ / c ( C18 ) 

in the first parentheses of Equation (C 17) is the fluctuating electrostatic potential in the 
frame of the parallel-drifting species s. 

Decomposing Sf± s into its adiabatic (Eq. C 17) and non-adiabatic (Ecp C 16) parts (cf. 
Antonsen & Lane 1980; Catto et al. 1981), the complete solution for f s may be written 
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Next we derive an evolution equation for the the gyrokinetic response h s . 


C.2.3. First Order, 0(e) 

Introducing the gyrokinetic potential 
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and using Equations (C 12) and (C 19), at this order Equation (C 1) becomes 
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(C 22) 


is the velocity of the guiding centre. Upon performing a ring average of Equation (C 21) 
over d at fixed R s , defined, for any function a(t,r, v), as 
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we find that the entire right-hand side of Equation (C21) vanishes. This follows from 
the periodicity of Sf 2 S in i9 and from the fact that, for any arbitrary function o(r), the 
ring average (t>x • V «)h s — 0 (see eq. A21 of Howes et al. 2006). Thus, the ring-averaged 
Equation (C21) is 
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Using the decomposition SB = VAm X z + SB\\z (Eq. C2) in Equation (C22) and 
retaining only first-order contributions, the ring-averaged guiding-centre velocity is 
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where we have used the identity (vx<5-B||).r 0 = — (Vx(^x • -Ax))h s - 

Substituting Equation (C25) into Equation (C24), we obtain the gyrokinetic equation 
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where the Poisson bracket is defined in the usual way: 
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The differences between Equation (C 26) and the usual (Maxwellian) slab gyrokinetic 
equation (cf. eq. 25 of Howes et al. 2006) lie entirely on the right-hand side, which arises 
from (ring-averaged) changes in the kinetic energy and magnetic moment of the parti¬ 
cles. In particular, the final term in Equation (C26), absent in Maxwellian gyrokinetics, 
ensures that h s evolves in such a way as to preserve adiabatic invariance. Indeed, if we 
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augment the gyrokinetic response in the following way,f 

h s = h s + ^^( x ) Ra , (C28) 

1 ±0s 

the gyrokinetic equation (C 26) can be expressed more compactly as 


dh s dh s c ~ qj$ a 
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These matters are discussed in detail in Section C.4. 
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C.3. Field Equations 

The equations governing the electromagnetic potentials are obtained by substituting 
Equation (C 19) into the leading-order expansions of the quasineutrality constraint (Eq. C 3) 
and Ampere’s law (Eq. C 4). This procedure requires computing parallel moments of the 
perpendicular-differentiated equilibrium distribution function: we denote 
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all of which equate to unity for a drifting bi-Maxwellian distribution (Eq. 2.24)4 
To first order, O(e), Equation (C3) becomes 
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Since the field variables ip and Ay are functions of the spatial variable r, the velocity 
integral of h s must be performed at constant location r of the charges rather than at 
constant guiding centre R s . This introduces a gyro-averaging operation dual to the ring 
average defined in Equation (C23): 


(h s (t,R s ,v\\,v±)) r = j) ddh s (t,r + V * * ,v\\,v ±^ . 
Equation (C31) may then be written as 

0 = H ( is 
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Likewise, the parallel and perpendicular components of Ampere’s law become, respec- 


f By Equation (C14), we can add any gyrophase-independent (at constant guiding centre 
R s ) function to the gyrokinetic response h s and still satisfy the zeroth-order kinetic equation. 

J If there are no interspecies parallel drifts and if fos is symmetric about Wy = tty 0s , then 
wliOsCi = 0. 
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is the pressure anisotropy of species s augmented by the parallel ram pressure from 
equilibrium parallel drifts. Upon integrating by parts with respect to the gyroangle, 
Equation (C35) can also be written as 


( B 0 6B I, \ 

Vx Vx : (6P± + / -^r 1 J = 0, (C 37) 

where 

SP± = ^2 / d 3 ^ m s (v±v±h s ) r (C38) 

s 

is the perpendicular pressure tensor. The perpendicular component of Ampere’s law is 
therefore a statement of perpendicular pressure balance for the compressive fluctuations. 

Together with the gyrokinetic equation (C 26), the field equations (C 33)-(C 35) provide 
a closed system that describes the evolution of a gyrokinetic plasma with non-Maxwellian 
f 0s and parallel species drifts. It remains to show that, in the kp i: m e /mi -C 1 limit, 
these equations reduce to those of KRMHD. 


C.4. Meaning o/d/ ls Boitz and h s : (ii||,i>x) vs. (s S yVs) Coordinates 

We pause the derivation here to offer a few comments regarding our choice of velocity- 
space coordinates, which will hopefully aid the reader’s grasp of the physical content 
encapsulated in the Boltzmann response (Eq. C 17) and the gyrokinetic equation (C26). 
Just as in the main text (see §4.2), we have opted for analytical convenience by choosing 
to work with vh and v± as our velocity variables. If, instead, our focus were to be on 
physical insight, then arguably better variables would be the total particle energy in the 
parallel-drifting frame, 


. l | , „i2 , 

£ s = £os +£is = ^rn s \v - u \\0s z \ + QsT s , 

and the gyrophase-dependent part of the first adiabatic invariant, 


Vs = VOs + Vis 


m s v\ 
2B 0 


9^ 

B 0 



(C 39) 


(C40) 


written out to first order in the fluctuation amplitudes (e.g. Kruskal 1958; Taylor 1967; 
Hastie et al. 1967; Catto et al. 1981; Parra 2013). The ring-averaged time derivatives of 
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both of these quantities, taken along a particle orbit, are first-order in e, namely 

^ ( d d \ 

( £ s)r, = q s ( Ql +M I|0 ^ ) <x)h. ~ O(eu}£ 0s ), 

= f (I + “< §;) {x)r ■ ~ 

and thus are the same order as the gyrokinetic equation (see §C.2.3). 

In Section C.2.2, we stated without proof that the Boltzmann response (Eq. C 17) 
arises from the evolution of /o s under the influence of the perturbed electromagnetic 
fields. Using 


(C41) 
(C 42) 


fl = -Tiio, 


dfos 

de, 


Os 


and 


fos = -T± Os 


f 1 dfos 
\Bq dflQg 


df 0s \ 

d£ 0s J 


(C43) 


to transform the (vu , ujJ-derivatives into (eo s , /ro s )-derivatives, the sum of the equilib¬ 
rium distribution function and the Boltzmann response may be written in the following 
suggestive form: 


fos( £ Os, HOs) “b ^/ls,Boltz — fos “b QsP, 
= fos + £ 1 


/ df 0s 


Bo 


ip . 


V W A \\ \ dfo 


d£ 0 s 

dfos , dfos 

~-b His ^- 

O£0s C'MOs 

— fos{ £ Os + £lsj HOs + His) = fos{ £ s, Hs) 


dnos 


(C 44) 


i.e., it is the electromagnetically perturbed distribution function fo s (£ s ,~p s ) Taylor ex¬ 
panded about /os(£os, Hos) and written out to first order in the change in particle energy 
and the (gyrophase-dependent) change in the magnetic moment. Thus, the Boltzmann 
response does not change the form of the equilibrium distribution function if the latter 
is written as a function of sufficiently precisely conserved particle invariants. 

Using Equation (C44), we can absorb the Boltzmann response into fo s and write the 
total distribution function of species s (Eq. C 19) as 


fs — fos( £ s, Hs) + h s (t, R s ,£ s , Hs) + ^/2s + ■ • ■ ■ 

The gyro-averaged Vlasov equation at O(eujfo s ) then becomes 

dh 


dh s 

~dt 


/ ri \ '' ‘ /— \ ^ f Os /— \ ^ f Os 

{Rs)r ‘ • m, - ~ {£s)r ^- { ^ )r ^~ 


(C45) 


(C46) 


Using Equations (C41) and (C42) for the ring-averaged rates of change of £ s and JL S , 
respectively, it then becomes clear that our gyrokinetic Equation (C 26) is simply the 
ring-averaged Vlasov equation, ( f s )R s = 0, written to lowest order in e. 

It is often convenient and physically intuitive to subtract from JL S a first-order gyrophase- 
independent piece, namely 

Hs = Hs- < jt(x)rs- ( C47 ) 

Bo 

In doing so, we gain an order in the ring-averaged conservation property, ( Hs)r , ~ 
O(e 2 uiH 0 s)i making Hs the asymptotically conserved adiabatic invariant. The additional 
contribution to ~p s i n Equation (C 47) also ensures that, at long wavelengths, Hs reduces 
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Ps = 


m s v 
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Similarly, we define the energy variable 

£ s = e s - qsiv'sW, 


v± ■ A a 


Rs 


(C48) 

(C49) 


which reduces in the long-wavelength limit to the energy variable defined by Equation 
(4.10): 
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(C 50) 


i.e., it is the kinetic energy of the particle as measured in the frame moving with the uj| 0s 
and E X B drifts. 

If we then write f s = fos{£si Ps) + Sf s , it is straightforward to show by using Equations 
(C28), (C44), (C47), and (C49) that the following expressions are equivalent to leading 
order: 


Sfs — Sf s /os(^Os; POs) fosi^s i Ps) 

dfos _ qp_ 
dsos B 0 
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(C 51a) 
(C 516) 

(C 51c) 

(C51d) 

(C51e) 


Equation (C 51) is the k±pi ~ 1 generalisation of the perturbed distribution function 
Sfs defined by Equation (4.6), which prominently features in the generalised free energy 
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(Eq. 5.1) of KRMHD; i.e., it is the perturbed distribution function if fo s is taken to be 
a function of (e s ,/x s ) instead of (um,uj_l). 

At long wavelengths, the difference between the Boltzmann response and its ring av¬ 
erage that features in Equation (C51c) is 


^/ls,Boltz (^/ls,Boltz).R s — 


q s v± X z 
T±os n s 

q s v ± X z 
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IIOs 


= -u± ■ 
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6B± 


<p(r) - 

V 
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ft 


Os 


v\\A\\(r) _ M yo^l|( r ) 
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dfos , . 


f II 

JO s 

M ||0s ^ dfos 
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dvw 


(C 52) 


where in the third line we have identified (cf. Eq. 2.25) 


u ± 

6B± 

yfAnpo 


X V ± ip(r) = z X Vj_$(r), 

Bo 

-7^2 X Vj_A||(r) = z X Vj.^(r). 
Bq 


(C 53a) 
(C 536) 


Comparing Equations (C51c) and (C51d), we see that the Alfvenic fluctuations com¬ 
prise the piece of the gyrokinetic response h a that is cancelled at long wavelengths by 
the Boltzmann response. Indeed, by substituting Equations (C52) and (C53) into the 
expression for the full distribution function, f s = fo s + <5/i Sj Boitz + h s + ..., we find that 
we may absorb the Alfvenic fluctuations into the equilibrium distribution: 


, / 1 | / -12 m s v± 

fo ° I 


fos ( \rn a {v\\ - u|| 0s ) 2 + \m s w \, 


m s w± 
2 Bn 


(C 54) 


In other words, Alfvenic fluctuations do not change the form of the distribution function, 
but rather define the moving frame in which any changes to it are to be measured. Phys¬ 
ically, this is because particles in a magnetised plasma adjust on a cyclotron timescale to 
take on the local E X B velocity. This principle is what underlies Kulsrud’s formulation 
of KMHD, in which the perpendicular particle velocities are measured relative to the 
E X B drift, the latter being governed by a set of MHD-like fluid equations rather than 
a kinetic equation. The implication—that discussed at length in the main text of this 
Paper (§2.5.3, in particular)—is that the compressive component of the turbulence is 
passively advected by the Alfvenic fluctuations in the inertial range. 


C.5. A Tactical Step: From Rings to Gyrocentres 

Considering the content of Equation (C 54)— that the Alfvenic fluctuations have a gyroki¬ 
netic response that is largely cancelled at long wavelengths by the Boltzmann response— 
it is often advantageous to not work directly with h s . This choice is made by most 
(Maxwellian) 8f gyrokinetic codes (e.g., Astro-GK ; Nurnata et al. 2010), which avoid the 
numerical error arising from this near-cancellation by working instead with the perturbed 
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distribution function 


q s / v ± ■ A ± 

g s = h s - — ( ip - 

Tos \ c 
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= n s + {ojis Boitz/R s + \ - 
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= (Sfu) Rs + ^-(^^) fos. 

To S \ C / Rs 


fos 


Rs 


In fact, many standard treatments of gyrokinetics use g s instead of h s . In the electrostatic 
limit, the use of g s (which, in this limit, equals (Sfi s ) Rs ) aids in the interpretation of 
polarisation effects within gyrokinetics (Krommes 2012), places the gyrokinetic equation 
in a numerically convenient characteristic form (Lee 1983), and arises naturally from the 
Hamiltonian formulation of gyrokinetics (Dubin et al. 1983; Brizard & Hahm 2007). 

We follow these practices and introduce, for non-Maxwellian /q s , 


q s / , tipAi, , 

g s = h s - —- ( ip s -) /, 
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II Os 


Os 


R s 


— h s + (<5/i s ,Boitz)ji s + — ( - ) fos, 

J- _LOs \ ^ 


R s 
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q s / v± • Aj 
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±0s 


fos, 


Rs 


(C 55a) 
(C 556) 
(C 55c) 


which is the k±p s ~ 1 generalisation of the perturbed distribution function g s defined 
in KRMHD (cf. Eq. 2.44). The physical distinction between the two formulations of 
gyrokinetics—one expressed in terms of h s , the other in terms of g s —is as follows. Writ¬ 
ten in terms of h Sl the set of gyrokinetic-Maxwell equations describes the dynamics of 
physically extended rings of charge as they move in a vacuum. If the fluctuating plasma is 
instead described by g Sl these equations describe a gas of point-particle-like gyrocentres 
moving in a polarizable medium (see Krommes 1993, 2012; Abel et al. 2013). Hence, we 
refer to g s as the gyrocentre distribution function. 

Using Equation (C55) to replace h B in the gyrokinetic equation (C26), we find that 
g s evolves according to 


dg s 

dt 


+”ii If + i, {(x)R - ' s - } = " “W { ' 0>R - ) 

) ft (O 56) 


1 9A\\ j. - / v± • A l 

c~dt +b ' Vip ~ b ’ V \ c 


Rs 


R s 


where 




d SB L 8 l r . 1 


Bo -- 82 Bo 

is the spatial derivative along the perturbed magnetic field (cf. Eq. 2.31).f As in S09, we 
have used compact notation in writing out the nonlinear terms: ({Am ,g> — (ip)r s })r s = 
({A||(t'),v?(t')})r s -{(A||) Rs , (ip)r s } 7 where the first Poisson bracket involves derivatives 
with respect to r and the second with respect to R s . 

The field equations (C33), (C34), and (C37) are best written in Fourier space. 


f When taking the gradient of a ring-averaged function, e.g. the term b- V(i>i • Ax)r s in 
Equation (C56), the velocity variable v is held fixed. 
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C. 6 . Fourier Space 

The ring-averaging (Eq. C 23) and gyro-averaging (Eq. C 32) procedures take on a rather 
compact form when expressed in Fourier space. First consider the gyrokinetic potential, 
decomposed into plane waves: \{t, r, v ) = Xk(t, v) exp(ifc • r). The ring average of its 
Fourier coefficient is 


(. Xk{t,v)) Rs = — $ M 


Fk 


w ll 


An 


V± A 


L k 


— do(Us) f (pk 


V\\ 


c 

All k 


exp 


—ifc 


v± X z 


Tj_os 2v]_ Ji(a s ) SB\\ k 


Vs <h_L s 


Bn 


(C 58) 


where a s = _/f2 s and we have used the definition SB = z • (ifc X Aj_*.). Sim¬ 

ply put, ring averaging amounts to multiplication by either the zeroth- (Jo) or first- 
order (Ji) Bessel function, depending on whether the additional angular dependence of 
v± appears in the integrand. Similarly, we can Fourier decompose g s (t, R s ,v\\,v±) = 
Efc 9sk{t, U||, tq_) exp(ifc • R s ) and perform its gyro-averages at constant r, viz. 


(gsk{t,v\\,v±)) r 
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27T 
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(C 59) 


(v±g s k(t,v\\,v±)} r = ^ <j> d^vj_g sfc exp ^i k- 
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(C 60) 


{v±v±g sk (t,v\\,v±)) r = ^ <j> ddt)j_t; ± g sfc exp ^ifc ■ 
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These operations transform integro-differential equations in real space into integro-algebraic 
equations in Fourier space. 


C.7. Field Equations 

In order to write the field equations in terms of g s concisely, we first define several vy-, 
v±_-, and Bessel-function-weighted integrals over the equilibrium distribution function, 
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suitably normalised: 
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(C 62a) 
(C 626) 

(C 62c) 

(C62 d) 
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1 — -a s Cn s + ... 

(C 62/) 
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where a s = k~ L p 2 s /2. In doing so, we have introduced two additional coefficients, given by 


Gn, = 


- 1 [d 3 vM^^Mf 0s and C 20s = — f d 3 v^M^f 0s , (C63) 
j v th±s u \\0s n 0s J 2 v? h±s 


^0 s 
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both of which equate to unity for a drifting bi-Maxwellian (Eq. 2.24).f Note the num¬ 
bering scheme used for the F^ m subscripts, which reflects the number of powers t of v 
and m of un in the integrand. 

With these definitions in hand, and working in the Fourier domain, we can express 
Equations (C33), (C34), and (C37), which represent respectively the quasineutrality 
constraint and the parallel and perpendicular components of Ampere’s law, as follows: 
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f If there are no interspecies parallel drifts and if fo s and / 0s are symmetric about wy = wji 0s , 
then then u'| 0 a Cu a = Wy 0 s r «i( a s) = 0 . 
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Note that ip k always enters the field equations in combination with —u'^ 0s A\\ k /c (cf. Eq. C 18). 


C.8. Massless Electron Fluid 

In this Section, we carry out an expansion of the electron gyrokinetic equation in powers of 
(we/biTi ) 1 / 2 ~ 0.02 (for hydrogen plasma). This expansion is done while still considering 
/3|| s , /3j_ s , T|| s , tj_ s , and k±pi to be order unity, so that k±p e ~ k±pi (m e /mi ) 1 ^ 2 <C 1. 
Then we can expand the Bessel functions arising from averaging over the electron ring 
motion, 

Jo{a e ) = 1 — ja 2 + ..., -= 1 — -i a 2 + ..., (C 67) 

4 a e 8 

and evaluate the fields <p, An, and SB y at r = R e . The electron kinetic equation, accurate 
up to and including the first order in (m e /m,y 1/2 , then reads 
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(C 68 ) 


where we have indicated underneath each term the lowest order to which that term enters 
when compared with vy dSf e /dz. We refer the reader to §4.1 and, in particular, equations 
(80)-(82) in S09 for details on obtaining this ordering. 

We expand g e = g + gi 1 ^ + ... in powers of (m e /mi) 1 / 2 and carry out the expansion 
to the first two orders. 
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C.8.1. Zeroth Order 

To zeroth order, the electron kinetic equation is 
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Equation (C69) is identical to the electron kinetic equation obtained from KRMHD 
(Eq. 2.36), as it should be (since a e -C 1). Multiplying Equation (C69) by T^ 0e /en 0e and 
integrating over the velocity space, we obtain 
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which matches the expression for (minus) the parallel electric held given by Equation 
(2.37). Substituting Equation (C70) back into Equation (C69) gives 
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from which follows the equations of state for the electrons, Equations (2.41) and (2.42). 

C.8.2. First Order 
At first order, Equation (C 68 ) reads 


„( 0 ) 


d9e +uyb- Vg A) + J- \ ip- 


dt 


T±0e vj_ SB\\ 


e «thJ_e Bo ,9e 


( 0 ) 


(C 72) 


21 


11 Oe 


l \\0e 


1 9 Ah j. Tj_ 0e t — v\ SBu 

c dt e <_Le Bo 


f H 
J 0e~ 


Using Equations (C 70) and (C71) and integrating Equation (C 68 ) over the velocity 
space, we find 
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where the parallel electron velocity is first order: 
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The hrst two terms in Equation (C 73) can be combined upon identifying 

d d dc 

as the Lagrangian time derivative measured in a frame transported at the E X B drift 
velocity, ue = —cV ]_T> X z/B 0 (cf. Eq. 2.30): 
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As we will show in the following section, Equation (C76) indicates that, at hrst order in 
the mass-ratio expansion, the magnetic hux is not tied to the E X B how. This marks a 
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departure from KRMHD, where k±pi is effectively zero at all relevant orders. However, 
Equation (C76) does not signal the breakdown of magnetic-flux conservation. 


C.8.3. Magnetic-Flux Conservation 

While it may not be readily apparent, Equations (C70) and (C76) state that the 
magnetic flux is frozen into the electron flow velocity and is therefore exactly conserved 
in that frame. To see that this is indeed the case, we begin by combining Equation 
(C 70) for the parallel electric field with the lowest-order expression for the perpendicular 
electric held — Vj jp' e associated with the Alfvenic fluctuations measured in the frame of 
the drifting electrons to obtain the total electric field, 
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being the effective velocity into which the magnetic flux is frozen. Therefore, there is a 
frame in which the magnetic flux is exactly conserved. We now show that this frame is 
associated with the electron flow velocity. 

We use Equation (C57) to rewrite Equation (C76) in the following equivalent form: 


8[ , d_ 
dt +U ^ 0e dz 


ABy 

B 0 


Sn e 

noe 


+ 


Bn 


Fe 


c\ e Tj_ 0e 8B\\ SB || 
B 0 


(7" e 

°0e 


Bo 


= b-V 


U\\ e 


(C80) 

To the left-hand side of this equation we add zero, written in a rather auspicious guise: 
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Equation (C 80) then becomes 
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and we may identify the first term in the Poisson bracket as the perpendicular component 
of the electron flow velocity, 
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It is clear that the leading-order contribution to M e ff is precisely that given by Equation 
(C 82). This identification made, we can interpret Equation (C 76) as the reduced electron 
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continuity equation. To wit, Equation (C81) expressed in terms of u ± e , 
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may be combined 
to O(eujB 0 ), 
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(C 83) 


with the parallel component of the induction equation (C 78) expanded 


( d d \ ^ jz> 

~Ql + M ||0e^7 + M -Le • ^-L J = & ' ^ u \\e ~ V • U e , (C 84) 

to find 

(s + "i'4 + “-- v ")^=- v '“- (C85) 

the reduced electron continuity equation. 


C.8.4. Field Equations 

The two fluid-like Equations (C 70) and (C 76) form the system that describes the 
electrons. These are closed by the gyrokinetic equation for gi (Eq. C 56 with s = i), and 
by the three integral relations derived from quasineutrality (Eq. C 64) and the parallel 
(Eq. C 65) and perpendicular (Eq. C 66) components of Ampere’s law. We now express 
the latter explicitly in terms of Sn e , u\\ e , ip, Ay, 6B\\, and c/j. 

Expanding the summation over species and using Equations (C 71) and (C 74) to com¬ 
pute the velocity-space moments of g e , Equations (C64), (C65), and (C66) become, 
respectively, 
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where we have used JU c i u '\\oi = u \\o e — 0 (which follows from the mass-ratio expansion). 
The lowest order in k±pi at which each term enters is indicated underneath that term, 
following the subsidiary ordering discussed in §5.2 of S09. As promised, if we retain 
only the zeroth-order terms, Equations (C 86)-(C 88) reduce to their respective KRMHD 
Equations (2.48)-(2.50). 

C.9. Inertial-Range Turbulence: Reduction to KRMHD 

Thus far, we have constructed the theory for electrons, which determines the equations 
of state of the electron fluid, evolves the parallel component of the vector potential via 
a generalised Ohm’s law, and demonstrates that the magnetic flux is convected by the 
perpendicular electron flow. In this Section, we proceed to derive the gyrokinetic theory 
for the ions, and show that it reproduces the KRMHD equations (2.51a)-(2.51c) in the 
long-wavelength limit relevant to the inertial range. 

C.9.1. Compressive Fluctuations 

Substituting the expression for dA\\/dt that follows from Equation (C70) into the 
gyrokinetic equation (Eq. C56), we find that gi satisfies 
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/ Ri 


We have indicated underneath each term the lowest order in k±pi at which that term 
enters. The zeroth-order terms should reduce to the ion drift-kinetic equation (Eq. 2.51c). 
Using Equations (C 75) and (C 57) to group terms on the left-hand side of Equation (C 89) 
into the Lagrangian operators d/d t and b ■ V, and replacing (Z i e/T\\ 0 i){v± • A±/c) Ri with 
the lowest-order expression for it, — {v‘]_/vl^ i )(SB\ t JB 0 ), we find 
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as promised. 
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C.9.2. Alfvenic Fluctuations 

If we now multiply Equation (C89) by Ci/not, sum over the ionic species, integrate 
over the velocity space (keeping r constant), and use Equations (C86) and (C87) to 
express the velocity-space integrals of gi, we obtain 
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It is straightforward to show that the zeroth-order component of this equation is identical 
to Equation (C 73) and so vanishes. Next we consider the first-order terms. Noting that 
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we multiply Equation (C 91) by —eno e /poj inverse Fourier transform back into real space, 
and use Equation (C 53) to get 
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The final term in parentheses vanishes by the reduced induction equation (2.27), leaving 
the reduced vorticity equation (2.29). 

This completes our derivation of KRMHD from the k±pi <C 1 limit of the non- 
Maxwellian gyrokinetic theory. 


Appendix D. Coefficients for a Bi-Kappa Distribution Function 

A bi-kappa distribution function is often used to describe the non-tlrermal electron 
population in the solar wind and, in particular, its suprathermal (T e ~ 60 eV) halo (e.g., 
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Vasyliunas 1968; Maksimovic et al. 1997a, b, 2005). In this Appendix, we evaluate all of 
the dimensionless C/ s coefficients introduced in the main text and in Appendix C for the 
bi-kappa distribution function 
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where T is the Gamma function, ac > 3/2 is the spectral index, and 
e \\s = 


and d ±s = v th±s \l - 


2 K 


(Dl) 


(D 2) 


are the effective parallel and perpendicular thermal speeds, respectively; u t h|| s and u t h± s 
are defined as in Equation (2.15). At low and thermal energies, the bi-kappa distribution 
approaches a Maxwellian distribution, whereas at high energies it exhibits a non-thermal 
tail that can be described as a decreasing power law. Note that Equation (D 1) tends to 
the bi-Maxwellian distribution (Eq. 2.24) as k — > oo. 

Evaluation of the C/ s coefficients is eased by rewriting /bi- K , s in integral form: 
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and from Equation (C63): 
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The coefficients defined by Equation (B 2) involve Landau-like integrals, which may 
be written in terms of the modified plasma dispersion function 
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where = (w — &|| w f|oi)/l^|| l^thU* is the dimensionless Doppler-shifted phase speed. 


Appendix E. Nomenclature 

In this Appendix, for the reader’s benefit we provide a glossary of frequently used 
symbols in our formulations of non-Maxwellian KRMHD and gyrokinetics. Each symbol 
is accompanied by a textual description and a numerical reference to either the section(s) 
in which the symbol was introduced or, if available, the equation by which the symbol was 
defined (given in parentheses). Throughout the manuscript, the subscript “0” appended 
to any of the following symbols denotes an equilibrium value; the pre-factor “<5” denotes 
a fluctuation. The species index s = i (for ion), e (for electron), or a (for alpha). 

Miscellaneous 


m s 

mass of species s 

§2.1 

Qs 

charge of species s (= Z s e) 

§2.1 

e 

electric charge magnitude 

§2.1 

c 

speed of light 

§2.1 

\ 

unit dyadic 

§2.1 

e 

expansion parameter (<C l) 

§2.3, §C.l 

Jo, Ji 

zeroth- and first-order Bessel functions 

§C.6 

Lengthscales 


L 

fiducial macroscale 

§ 1 , §c 

A m fp 

collisional mean free path (= u t h|| i/va) 

§1 

Ps 

gyroradius of species s (= "CthJ .s/^s) 

§2.1, §C 

k~ 1 

K\\ 

(inverse) parallel wavenumber 

§2.3, §C.l 

h 1 

h-‘ 

(inverse) perpendicular wavenumber 

§2.3, §C.l 

Frequencies 


OJ 

frequency of the fluctuations (= co T + iy) 

§2.3, §C.l 

Vii 

ion-ion collision frequency 

§ 1 , §c.i 

is 

dimensionless phase speed (= w/fc||i' t h||s — ^yos/^thlls) 

§4.4, §B.l 

12s 

gyrofrequency of species s (= q s Bo/m s c) 

§ 1 , §c 


Phase-Space Coordinates 
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v \\ 

velocity-space coordinate parallel to the magnetic field 

§2.1, (C 8) 

v± 

velocity-space coordinate perpendicular to the magnetic 
field 

§2.1, (C 8) 

w_ L 

velocity-space coordinate perpendicular to the magnetic 
field and peculiar to the mean perpendicular flow of 
species s 

§2.1 

r 

real-space coordinate 

§2.1 

£ s 

kinetic energy of a particle of species s as measured in 
the frame of the Alfvenic fluctuations and the equilibrium 
species drift 

(4.10), (C49) 

Ms 

first adiabatic invariant of a particle of species s 

(4.11), (C47) 

i9 

gyrophase angle 

(C 8) 

R s 

guiding-centre position of species s 

(CIS) 

e s 

total energy of a particle of species s in the frame of the 
equilibrium species drift 

(C 39) 

£o s 

kinetic energy of a particle of species s in the frame of 
the equilibrium species drift 

(C 39) 

& Is 

first-order correction to the kinetic energy of species s 

(C 39) 

Ms 

gyrophase-dependent part of the first adiabatic invariant 
of a particle of species s 

(C40) 

MOs 

lowest-order magnetic moment of a particle of species s 
(= m s vj_/2B 0 ) 

(C 40) 

Ml s 

first-order gyrophase-dependent correction to the mag¬ 
netic moment 

(C 40) 

Qj S 

dimensionless velocity-space coordinate perpendicular to 
the magnetic field (= k±v±/fl s ) 

§C.6 

Distribution Functions 


fs 

distribution function of species s 

§2.1 

f\\ 

J Os 

dimensionless derivative of the equilibrium distribution 
of species s with respect to the square of the parallel 
velocity peculiar to the equilibrium species drift 

(2.22), (C11) 

fos 

dimensionless derivative of the equilibrium distribution 
of species s with respect to the square of the perpendic¬ 
ular velocity 

(2.22), (C 11) 

fbi-M,s 

bi-Maxwellian distribution function of species s 

(2.24) 

9s 

perturbed distribution function if f s is taken to be a func¬ 
tion of uy and fi a (in KRMHD); gyrocentre distribution 
function (in gyrokinetics) 

(2.44), (C55) 


fs 

9 * 


h s 

3 fls,Boltz 

h s 


fhi-n,s 
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distribution function of species s as a function of (e s , p s ) 

passively mixed, undamped, ballistic component of the 
perturbed ion distribution function gi 

gyrokinetic response 

leading-order Boltzmann response 

gyrokinetic response corrected for g s conservation 

bi-kappa distribution function of species s 


Moments of the Zeroth-Order Distribution Function 

C'f s dimensionless coefficients related to perpendicular mo¬ 

ments of the parallel-differentiated equilibrium distribu¬ 
tion function; includes Landau resonance for ions 

Zy[ Maxwellian plasma dispersion function 

Cf- S dimensionless coefficients related to parallel moments of 

the perpendicular-differentiated equilibrium distribution 
function 

Ffm several dimensionless Vjf-, vTT-, and Bessel-function- 

weighted integrals over the (perpendicular-differentiated) 
equilibrium distribution function 

Z K kappa plasma dispersion function 


Densities 

n s number density of species s (= / d 3 v f s ) 

p volume density of plasma (= Y s m s n s ) 

c s charge-weighted ratio of number densities (= Z s no s /no e ) 

Velocities 

u s mean velocity of species s (= nj 1 J d 3 u u/ s ) 

m_l s mean perpendicular velocity of species 

u centre-of-mass velocity (= Y s m s n s u s / m s n s ) 

u± perpendicular centre-of-mass velocity (= cE X B / B 2 ) 

uj| s mean parallel velocity of species s measured in a frame 

comoving with the centre-of-mass velocity 

i>a Alfven speed 

Uth||s parallel thermal speed of species s 

I'thJ.s perpendicular thermal speed of species s 

va* effective Alfven speed 

d|| s effective parallel thermal speed of species s for a bi-kappa 

distribution function 


§4.2 

(4.50) 

(C 16) 

(C 17) 

(C 28) 

(Dl) 

(2.38), (B2) 

(4.26) 

(C 30) 

(C 62) 

(D6) 

§ 2.1 

§ 2.2 

§2.3 

§ 2.1 

§ 2.1 

§ 2.2 

§ 2.2 

(2.9) 

(2-14) 

(2.15) 

(2.15) 

(2.34) 

(D2) 
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9± s effective perpendicular thermal speed of species s for a (D2) 

bi-kappa distribution function 

Pressures 


Ps 

pressure tensor of species s 

(2.4) 

P\\s 

parallel pressure of species s 

(2.5) 

P±s 

perpendicular pressure of species s 

(2.6) 

P\\ 

parallel pressure of plasma (= E s P||s) 

§2.2 

Pa l 

perpendicular pressure of plasma (= Es^-Ls) 

§2.2 

A s 

dimensionless pressure anisotropy of species s 

(2.20) 

/3|| s 

ratio of parallel pressure of species s to the magnetic 
pressure 

(2.21) 

P±s 

ratio of perpendicular pressure of species s to the mag¬ 
netic pressure 

(2.21) 

Al 

ratio of parallel pressure of plasma to the magnetic pres¬ 
sure (= E S /3|| S ) 

§2.3 

P± 

ratio of perpendicular pressure of plasma to the magnetic 
pressure (= E s P±s) 

§2.3 


dimensionless pressure anisotropy of the electrons 
weighted by C\ a 

(2.39), (B4) 

A s 

dimensionless pressure anisotropy of species s augmented 
by the parallel ram pressure from equilibrium parallel 
drifts 

(C 36) 

Temperatures 


T\\s 

parallel temperature of species s 

(2.5) 

T ±s 

perpendicular temperature of species s 

(2.6) 

T ||s 

ratio of parallel temperature of species s to the parallel 
electron temperature 

(2.19) 

T±s 

ratio of perpendicular temperature of species s to the 
perpendicular electron temperature 

(2.19) 

Electromagnetic Fields and Potentials 


E 

electric field 

§2.1 

B 

magnetic field 

§2.1 

b 

unit vector in the magnetic-field direction (= B/B) 

§2.1 

E\\ 

parallel electric field 

(2.8), (2.37) 

E± 

perpendicular electric field (= — u± X B/c) 

§2.2 

<f> 

velocity stream function 

(2.25) 

4/ 

magnetic flux function 

(2.25) 
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C ± 

v 

A 

3 

P's 


X 


generalised Elsasser potentials (3.3) 

generalised Elsasser fields (3-4) 

electrostatic scalar potential (C 2) 

magnetic vector potential (C 2) 

current density (C4) 

electrostatic scalar potential in the frame of the parallel- (C 18) 

drifting species s 

gyrokinetic potential (C 20) 


Functions Defined for Parallel Kinetics (§4.3) 


G n 

integral over the perpendicular velocity space in 
Eq. (2.51 d), which relates the fluctuating electron num¬ 
ber density to the lower-order moments of g t 

G b 

integral over the perpendicular velocity space in 
Eq. (2.51e), which relates the fluctuating magnetic-field 
strength to the lower-order moments of gi 

4 

various perpendicular moments of the parallel- 
differentiated equilibrium ion distribution function 

Ai~> 

coupling coefficients quantifying the influence of the den¬ 
sity and magnetic-field-strength fluctuations on the ki¬ 
netic fluctuations; superscripts are nn, nB , Bn, and BB 

G ± 

eigenvectors resulting from diagonalising the reduced ion 
kinetic equation for a bi-Maxwellian plasma 

A ± 

(inverse) eigenvalues resulting from diagonalising the re¬ 
duced ion kinetic equation for a bi-Maxwellian plasma 

<?i 

useful ion coefficient 

Q 

useful ion coefficient 

TUi 

useful ion coefficient 

«< 

useful ion coefficient 


(4.15) 

(4.15) 

(4.17) 

(4.18) , §A 

(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 
(4.42) 


Invariants 

If A\v Alfven-wave invariants associated with ^ 

Waw Alfven-wave invariant 

Wcompr compressive invariant 

Wcompr compressive invariants associated with G ± 

Wg t compressive invariant associated with 7j, 

W generalised free energy 


(3.9) 

(3.10) 
(4.4) 
(4.37) 
(4.52) 
(5.1) 


56 M. W. Kunz, A. A. Schekochihin, C. H. K. Chen, I. G. Abel, & S. C. Cowley 


Differential and Integral Operators 


D/D t 

time derivative measured in a frame co-moving with the 
Alfvenic fluctuations and streaming along the magnetic 
held at velocity dm (= d/dt + u± s • V -1- v\\b • V) 

§2.1 

£ 

differential operator measuring the velocity-space 
anisotropy of a distribution function 

(2.23), (C 11) 

{•, ‘I 

Poisson bracket 

(2.28), (C27) 

d/di 

time derivative measured in a frame co-moving with the 
Alfvenic fluctuations 

(2.30), (C75) 

b • V 

space derivative measured along the exact magnetic-field 
direction 

(2.31), (C57) 

fa 

time derivative of the function f s (t, r, v) taken along the 
full phase-space trajectory of a particle of species s 

(Cl) 


ring average over i? at fixed R s 

(C 23) 

<--)r 

gyro-average over i9 at fixed r 

(C 32) 
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